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We propose a new discrete model — the twisted quantum double model — of 2D topological phases 
based on a finite group G and a 3-cocycle a over G. The detailed properties of the ground states 
are studied, and we find that the ground-state subspace can be characterized in terms of the twisted 
quantum double D a (G) of G. When a is the trivial 3-cocycle, the model becomes Kitaev's quantum 
double model based on the finite group G, in which the elementary excitations are known to be 
classified by the quantum double D(G) of G. Our model can be viewed as a Hamiltonian extension 
of the Dijkgraaf-Witten topological gauge theories to the discrete graph case with gauge group 
being a finite group. We also demonstrate a duality between a large class of Levin- Wen string-net 
models and certain twisted quantum double models, by mapping the string-net 6j symbols to the 
corresponding 3-cocycles. 
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I. INTRODUCTION 

The study of possible phases of matter has gone beyond 
Landau's paradigm of symmetry breaking for decades, 
which leads to the discovery of topological phases of mat- 
ter. Among all possible topological phases, there are a 
class of them that are believed to bear intrinsic topologi- 
cal order—, in which they display features such as robust 
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ground state degeneracies, Abelian or non-Abclian braid- 
ing (anyonic) statistics of quasi-particle excitations, and 
in many cases protected edge excitations. The classic 
examples of these phases include the (fractional) quan- 
tum Hall states, Z2 spin liquids, chiral spin liquids, and 
p + ip superconductors^—. The physical characteristics 
of topological phases urge the search of the mathematical 
structures that classify the topological phases. It is then 
natural to resort to certain theoretical models that can 
yield various topological phases. 

There is a very general framework — the string-net 
models^, also known as the Levin-Wen models — 
supplying exactly soluble models that incorporates a 
large class of intrinsically topological phases, notably 
those preserving time-reversal symmetry. Although it is 
believed that tensor category theory is the mathematical 
framework that underlies these models, a general classi- 
fication of these models — in particular of the topological 
phases they describe — is yet to be found. 

The intrinsically, topologically ordered systems are 
roughly speaking those gapped quantum phases of mat- 
ter that involve long range entanglement (LRE). In con- 
trast, there are gapped quantum phases of matter that in- 
volve short range entanglement (SRE), which, when sym- 
metry is unbroken, give rise to nontrivial phases, called 
symmetry-protected topological (SPT) phase o 12 ' 13 , such 
as the Haldane phase on one-dimensional spin chainis 
and topological insulators^—. Characteristic proper- 
ties of these phases are usually non-degenerate ground 
states and, if the system has a boundary, nontrivial edge 
excitations. 

Very recently, however, it is discovered that a specific 
SPT phase, namely an Ising spin model with a gauged 
Z2 symmetry, admits a dual LRE phase described by 
a string net model^i. This remarkable duality is then 
conjectured^ to exist between a general SPT phase with 
discrete, gauged symmetry G and a string net model 
with fusion rules also given by the product rule of G. 
Soon after, this conjecture is confirmed in Re^, which 
henceforth implies that the classification of SPT phases 
provided by group cohomology in 2+1 dimensions via 
H 3 {G,U{1)) described in Re£&2i 43 m di rec tly provide 
classifications of the corresponding string net models. 

This classification of string-net models seems feasible, 
as the building blocks of these models, namely the 6j- 
symbols may fall into equivalence classes that are related 
to the 3-cocycles in the cohomology group H 3 (G,U(1)) 
of the symmetry group G of the model^. Nevertheless, 
in the string-net models that have been studied so far, 
the 6j symbols are assumed to respect the full tctrahc- 
dral symmetry, which may be too restrictive for a de- 
scription of topological phases. Namely, as pointed out 
in ~Re&, the topological phases described by the Levin- 



Re£ss offers a field theoretic approach that obtain the same clas- 
sification. 



Wen model with tetrahcdral symmetry may not account 
for all topological phases classified by H 3 {G, U(l j). 

This has motivated us to propose a new class of dis- 
crete models for 2D topological phases, called the twisted 
quantum double model for reasons to be clear later, 
whose construction involves a 3-cocycle, an element in 
the group cohomology group H 3 (G,U(1)). More pre- 
cisely we consider a model on a planar graph of triangles, 
each edge of which is graced with a group element of a 
finite group G. The Hamiltonian of the model has ma- 
trix elements constructed with a 3-cocycle a belonging 
to the cohomology group H 3 (G, U(l)) of G. We require 
that a satisfies only the 3-cocycle condition 8a = 1 where 
6 is the coboundary operator, which under the circum- 
stance of this paper is actually the pentagon identity in 
disguise. Owing to the absence of extra conditions put 
in by hand on a, All solutions to the 3-cocycle condition 
but one — namely the trivial 3-cocycle — do not respect 
the tetrahedral symmetry. In other words, any element 
of H 3 (G, U(l)) defines an instance of our model. 

We study our model in detail by placing it on a torus. 
In terms of 3-cocycles we construct, for the ground states 
of our model, explicitly three topological obscrvablcs, 
namely, the ground state degeneracy (GSD), the S and 
T operators that are a representation of the generators of 
the modular group SL(2, Z). This construction is a new 
result of ours, purely based on our model and in terms 
of the 3-cocycles of G, without using group representa- 
tion theory. These topological observables on the ground 
states lead to a set of topological numbers, respectively 
formed by the GSD, elements of S- matrix, and topologi- 
cal spin. We show that these topological numbers depend 
on the cohomology classes [a] e H 3 (G, U(l)). Moreover, 
equivalent 3-cocycles define equivalent twisted quantum 
double models, in the sense that their Hamiltonians can 
be continuously deformed into each other. Unfortunately, 
in this paper, we are not able to offer an answer to 
whether two inequivalent 3-cocycles always yield two dis- 
tinct topological phases. Instead, we present a few char- 
acteristic properties of the topological numbers, which 
may help to resolve this open question in future work. 
On top of these abstract constructions, we work out a 
few concrete examples for certain finite groups, Abelian 
and non-Abelian. 

We also discourse on how our model relates to topo- 
logical field theories and models of topological phases. 
It turns out that our model is a reasonable Hamiltonian 
extension of the Dijkgraaf-Witten theory2£r— of topo- 
logical Chern-Simons gauge theory in three dimensions, 
as we can identify the ground states of our model de- 
fined by an [a] s H 3 (G,U(1)) on the boundary of a 
three-manifold with the gauge-invariant boundary states 
of the Dijkgraaf-Witten theory defined by the same 
[a] in the bulk, which then equates the GSD of our 
model with the partition function of the correspond- 
ing Dijkgraaf-Witten theory. Since three-dimensional 
topological Chern-Simons theory corresponds to two- 
dimensional rational conformal field theory (RCFT)2£, a 
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connection between our model and RCFT is thus estab- 
lished. In particular, the GSD of our model with group 
G agree with the number of primary fields in the RCFTs 
that are orbifolds by the symmetry group G of a holo- 
morphic CFT. 

We demonstrate that our twisted quantum double 
model reduces precisely to Kitaev's quantum double 
model in the special case where the defining 3-cocycle 
is trivial. The nontrivial 3-cocycles in our model may 
twist the usual group algebra C[G] into a twisted group 
algebra, which mainly motivates the name of our model. 

As our model is motivated by the Levin-Wen model, 
we demonstrate a duality between a large class of Levin- 
Wen string-net models and certain twisted quantum dou- 
ble models, by mapping the string-net 6j symbols to the 
corresponding 3-cocyclcs. 

We would like to insert as an aside here that we min- 
imized the complexity of the mathematics in this paper 
without sacrificing the preciseness and comprchcnsibility 
of our presentation. For instance, although group coho- 
mology is a key word of this paper, but we assume zero 
prior knowledge of it because we define and present the 
n-cocycles as merely U(l) functions that satisfy an alge- 
braic condition. As such, we believe the paper is acces- 
sible to a wide range of physicists and mathematicians. 

Our paper is organized as follows. In Section |TT] we 
construct our new model of topological phases. Section 
TTll is devoted to the general setting for the topological 
observables. In Section HVl wc compute the ground state 
degeneracy (GSD) on a torus and study the correspond- 
ing topological degrees of freedom. Section |V] furnishes 
the construction of two more topological observables that 
give rise to fractional topological numbers. We present 
a classification of the topological numbers in our model 
in In Section IvTl Section IVlII offers concrete examples of 
our model for a number of finite groups. The next three 
Sections I Villi HXl and IXl relate our model respectively to 
Kitaev's quantum double model, Dijkgraaf-Wittcn topo- 
logical gauge theory, and Levin-Wen string-net model. 
The final section (Section IXIjl concludes with remarks 
and outlook. 



II. THE MODEL 

In this section, we shall construct our model in (2 + 1)- 
dimension, as an exactly-soluble Hamiltonian on the 
Hilbert space spanned by planar graphs consisting of tri- 
angles whose edges are graced with group elements in 
certain finite group. 



A. Basic Ingredients 

The model is defined on a two-dimensional graph Y 
consisting of triangles only (Fig. [T]). Such a graph does 
not have any open edge and may be thought as a simpli- 
cial triangulation of certain two-dimensional Ricmannian 



surface, e.g., a sphere; however, in this model, we shall 
take the graph as abstract without referring to its topo- 
logical background except when we compare the model 
with other models, such as Dijkgraaf-Witten discrete 
topological gauge theories. Note that Fig. Q] is a crop 
of one such graph, so the open edges in the figure are not 
really open. We enumerate the vertices of T by any or- 
dered set of labels. The enumerations of the vertices we 
choose is irrelevant as long as their relative order remains 
consistent during the calculation. 




FIG. 1: A portion of a graph that represent the basis vec- 
tors in the Hilbert space. Each edge carries an arrow and is 
assigned a group element denoted by [ab] with a < b. 

The model is characterized by a triple (H, G, a), which 
can be denoted by Hc,a f° r short. The first in the triple 
is the Hamiltonian H. The second ingredient G is a finite 
group. Each edge of T is graced with a group element of 
G. The Hilbert space is spanned by the configurations of 
group elements on the edges of T. Each edge (see Fig. [1]) 
carries an arrow that goes from the vertex with a larger 
label to the one with a smaller label. To each edge e of 
the graph T, we assign a group element g e € G, and all 
possible assignments form the basis vectors of the Hilbert 
space. 



{5i,52, -,9e} 
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where E is the total number of edges in T. 

It convenient to denote both an edge and the group el- 
ement on the edge by simply [ab] with a < b the two 
boundary vertices of the edge. It is understood that 
[ba] = [ab]^ 1 . The inner product of the Hilbert space 
is the obvious one: 



=8[ab] [a'b'] 8 [be] [b'c']8[ac] [a'c'J 



(2) 

where only one triangle in T is drawn, and the ". . . " 
omits the (5-functions on all other triangles that are not 
shown. Note that three group elements on the three sides 
of a triangle, e.g., the [ab], [be] and [ac] on the RHS of 
Eq. ©, are independent of each other in general, i.e., 
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[ab] ■ [be] 7^ [ac] . From now on, we shall neglect the group 
elements on the edges but keep only the vertex labels 
when we draw a basis vector. 





FIG. 2: (a) The defining graph of the 3-cocycle 
a([viv 2 ], [v 2 v 3 ], [v 3 v 4 ]). (b) For a([viv 2 ], [v 2 v 3 ], [«3«4]) _1 . 

The third element is a normalized 3-cocycle a G 
H 3 (G, U(l)), i.e., a function a : G 3 -> 17(1) that sat- 
isfies the 3-cocycle condition 

a(gi,g2,93)a(g ■ 5i,52,33) _1 x (3) 
a(ffo,5i ' 52,53)«(5o,5i7ff2 ■ ff3) _1 a(5o, 9i, 92) = 1 
for all gi G G, and satisfies the normalization condition 



oiO-,g,h) = a(g,l,h) = a(g,h,l) = 1, 



(4) 



whenever <?, h G G arc arbitrary. We emphasize that this 
normalization condition is not an ad hoc condition we 
imposed as an extra on the 3-cocycles; rather, it is a 
natural condition that any group 3-cocycle can satisfy 
for the following reason. A 3-cocycle a is in fact an 
equivalence class of the 3-cocycles that can be scaled 
into each other by merely a 3-coboundary 6/3, where /3 
is a 2-cochain. It can be shown that for any equivalence 
class of 3-cocycles, there always exists a representative 
that meets the normalization condition in Eq. which 
is in turn justified. 

Note that every group has a trivial 3-cocycle ao = 1 on 
the entire G. One can define a 3-cocycle on any subgraph 
composed of three triangles, which share a vertex and 
any two of which share an edge. Consider Fig. 2(a) 
as an example: The four vertices are in the order V\ < 
V2 < V3 < V4; we define the 3-cocycle for this subgraph 
by taking its three variables from left to right to be the 
three group elements, [^1^2], [^2^3] and [^3^4], which are 
along the path from the least vertex v\ to the greatest 
vertex V4 passing V2 and v 3 in order; hence, the 3-cocycle 
reads a(k[i£2], [^2^3], [^3^4]). If one lifts the vertex V2 
in Fig. 2(a) above the paper plane, the three triangles 
turns out to be on the surface of a tetrahedron. In this 
sense, one can think of the 3-cocycle as associated with 
a tetrahedron as well, which is useful when the graph is 
really interpreted as the triangulation of a Ricmannian 
surface. 

On the other hand, if one switches the vertices V2 and 
V3 in Fig. 2(a)| one obtains Fig. 2(b)| which defines the 
inverse 3-cocycle ccQuit^]) [^2^3], [v3i>4]) _1 . Whether a 
graph defines a 3-cocycle a or the inverse a -1 depends 
on the orientation of the four vertices in the graph by the 
following rule. One first reads off a list of the three ver- 
tices counter-clockwise from any of the three triangles of 



2(b) If the list can be turned into ascending or- 
by even permutations, such as {vi,V2,V3,Vi) from 
2(a)| one has an a but an a" 1 otherwise, as by 



the defining graph of the 3 cocycle, e.g., (v 2 ,v 3 , V4) from 
Fig. 2(a) and (^3,^2, ^4) from Fig. |2(b)[ One then ap- 
pend the remaining vertex to the beginning of the list, 
e.g., (vi,V2 ,V3,V4,) horn Fig. |2(a)| and (vi,v 3 , v 2 , va) from 
Fig. " 
dcr 

Fig. 2(a)| one has an a but an a 
(vi,v 3 , v 2 ,Vi) from Fig. |2(b)| 

We would like to warn the reader of some abuse of lan- 
guage in the rest of the paper. For example, when we 
say "a 3-cocycle", we may refer to a class [a], a repre- 
sentative a, or the evaluation of a on a tetrahedron. For 
another example, although there is abstractly only one 3- 
cocycle condition as in Eq. ([3]), we may sometimes mean 
3-cocycle conditions by the evaluation of the condition 
on different tetrahedra. But all such and such should not 
cause any confusion contextually. 



B. The Hamiltonian 

The 3-cocycles will appear in the matrix elements of 
the model's Hamilton defined as follows. 



(5) 



where Bf is the face operator defined at each triangular 
face /, and A v is the vertex operator defined on each 
vertex v. We now elaborate more on these operators. 
The action of Bf on a basis vector is 



B, 




Ull'2j-Lf2l>3Hl>3l'lJ 




The discrete delta function 5\ 



VlV 2 ]-[v 2 V 3 ]-[v3Vl] 



(G) 

is unity if 



[^1^2] • [^2^3] • [^s^i] = 1, where 1 is the identity clement 
in G, and otherwise. Note again that here, the 
ordering of v\ , v 2 , and V3 does not matter because of 
the identities filv 1 v 2 }-[v2v 3 }-[v3v 1 ] = ^Iv3v 1 ]-[v 1 v2]-[v2v 3 ] 

and < 5['UiU2]-['U2'U3]-['!'3'Ul] < ^{['Ulf2]-[' u 2'l'3]-['l'3'l'l]} _1 

^WiI-M^sI-Mi-i^]- 1 = S [viv 3 ]-[v 3 v 2 ]-[v 2Vl ]- In other 
words, in any state on which Bf = 1 on a triangular 
face /, the three group degrees of freedom around v is 
related by a chain rule: 



[VlV 3 \ = [VlV 2 \ ■ [v 2 v 3 \ 



(7) 



for any enumeration v\ , v 2 , v 3 of the three vertices of the 
face /. 

The operator A v is a summation 



^ v ~ \Q\ E 

1 1 [vv']=geG 



(8) 



which deserves explanation. The value \G\ is the order 
of the group G. The operator acts on a vertex v with 
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a group element g € G by replacing v by a new enu- 
meration v' that is less than v but greater than all the 
enumerations that are less than v in the original set of 
enumerations before the action of the operator, such that 
[v'v] = g. A 9 , does not affect any vertex other than v but 
introduces a U(l) phase, composed of 3-cocycles deter- 
mined by v' and all the vertices adjacent to v before the 
action, to the resulted state. In a dynamical language, 
v' is understood as on the next "time" slice, and there 
is an edge [v'v] E G in the (2 + 1) dimensional "space- 
time" picture. Consider a trivalent vertex as an example 
(see Eq. ©). Without loss of generality, we assume that 
the enumerations of the four vertices are in the order 
v\ < v 2 < v 3 < V4. The basis vector on the LHS of Q 
is specified by six group elements, [^1^3], [1*2^3], [^3^4], 
[W1U4], [w2«i], and [^2^4]. The action of A 9 3 on this state 
reads 



v 3 



-S[v' 3 v 3 ], g a ([v!V 2 ], [v 2 v' 3 ], [v' 3 v 3 ]) a ([v 2 v' 3 ], [v' 3 v 3 ], [v 3 v 4 ]) 




x a([viv' 3 ], [v' 3 v 3 ], [v 3 v 4 ]) 1 




(9) 



where on the RHS, the new enumerations are in the order 
vi < v 2 < v 3 < v 3 < v 4 , and the following chain rule of 
group elements on the edges holds. 



[v' 3 vi] = [v' 3 v 3 ] ■ [u 3 ui], 
[v' 3 v 2 ] = [v' 3 v 3 ] ■ [v 3 v 2 ], 
[ v ' 3 v 4 ] = [ v ' 3 v 3 ] ■ [v 3 v 4 ]. 



The phase factor consisting of three 3-cocycles on the 
RHS of Eq. §§§ encodes the non-vanishing matrix ele- 
ments of Byl , namely 

V^JWJWMvt] M> 

=a ([v!V 2 ], [v 2 v' 3 ], [v' 3 v 3 ])a{[v 2 v' 3 ], [v' 3 v 3 ], [v 3 v 4 ]) (H) 
x a([viv' 3 ], [v' 3 v 3 ], [v 3 v 4 ]y ■ 

For each vertex on the LHS of Eq. Q, we group 
its three neighboring enumerations together with the 
new enumeration v' 3 in the ascending order. Hence, we 
have {vi,v 2 ,v' 3 ,v 3 ) for the lower vertex, (vi,v' 3 ,v 3 ,v 4 ) 
for the upper left vertex, and (v 2l v' 3 , v 3l v 4 ) for the 
upper right one, and then assign three 3-cocycles re- 
spectively to the three vertices: a ([^1^2], [^2^3], [^3^3]), 
a {[v 2 v' 3 ], [v' 3 v 3 ], [v 3 v 4 ]), and a {[v!V 3 ], [v' 3 v 3 ], [v^y 1 . 
The chirality of a 3-cocyle, or in other words, whether 
a vertex contributes a 3-cocycle a or the inverse a -1 , is 
based on the following criteria. We write down a triple 
for the three neighboring enumerations around each ver- 
tex in the counterclockwise direction and append v' 3 to 



the front, namely, (v' 3 , v%, v 2 , v 3 ) for the lower vertex, 
(v' 3 ,vi,v 3 ,v 4 ) for the upper left one, and {v' 3 ,v 2 ,v 4 ,v 3 ) 
for the upper right one. If it takes (odd) even number of 
steps to permute a list to the ascending order, the vertex 
contributes (the inverse of) the corresponding 3-cocycle 
in the action. 

The matrix elements in Eq. can be better mo- 

tivated and understood in the following way. One may 
think that the graph evolves in "time" under the driven of 
the Hamiltonian. Focusing on the vertex operator only 
and considering the A^ 3 in Eq. the action of the 

operator creates a new "time" slice by replacing the orig- 
inal vertex v 3 by v' 3 and connects the two vertices in 
the "time" direction. This scenario is shown in Fig. [31 
which is made three-dimensional (2 + 1) to illustrate the 
"spacetime" picture and the relation between our model 
and Dijkgraaf-Witten discrete topological gauge theory 
to be addressed in Section HXl 




FIG. 3: The topology of the action of A 3 V3 . 



As in Fig. [3J we can view the original three triangles 
on the LHS of Eq. Q asa tetrahedron v\v 2 v 3 v 4 and the 
three new triangles as another tetrahedron v\v 2 v' 3 v 4 , of 

(10) 

which the vertex v' 3 lies inside v\v 2 v 3 v 4 because of the 
ordering v' 3 < v 3 . Since v' 3 and v 3 are connected, there 
are three more tetrahedra generated effectively by the ac- 
tion of the vertex operator, namely v\v 2 v' 3 V3, v 2 v' 3 V3V 4 , 
and viv' 3 v 3 v 4 . It looks like that the original tetrahedron 
is split into four tetrahedra. This splitting of tetrahe- 
dron implies the three chain rules in Eq. (fTUI) . which 
then enables us to endow the three tetrahedra viv 2 v' 3 v 3 , 
v 2 v' 3 v 3 v 4 , and viv' 3 v 3 v 4 respectively with the three 3- 
cocycles a([v 1 v 2 ],[v 2 v 3 },[ v 3 v 3]), « ([ v 2v' 3 ], [v' 3 v s ], [v s v 4 ]), 
and a ([viv' 3 ], [v 3 v 3 ], [i>3V 4 ]) , following the rule shown 
in Fig. 02 

The operator A® 3 in Eq. ^ is just an identity operator 
if [W3W3] = 1, i.e., the identity in G. In fact, according to 
Eq. we have the following matrix element 



a ([viv 2 ], [v 2 v' 3 ], 1) a ([v 2 v 3 ], 1, [v 3 v 4 ]) 
x a ([viv' 3 ], 1, [W3W4]) _1 , 



(12) 



which is unity, by the normalization condition 

The vertex operator in Eq. can naturally extends 
its definition from a trivalent vertex to a vertex of any 
valence higher than three. The number of 3-cocyles in 
the phase factor brought by the action of A 9 V on a vertex is 
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equal to the valence of the vertex. The chirality of each 3- 
cocycle in the phase factor follows the criteria described 
in the previous paragraph. It is clear that A 9=1 = I by 
the discussion above. 

It can be shown that all Bf and A v are projection 
operators and commute with each other (see Appendix 
A). As a result, the ground states and all elementary 
excitations arc thus simultaneous eigenvectors of all these 
local operators. Moreover, the elementary excitations are 
identified as local quasi-particles that are classified by 
the the representations of the local operators. 

We shall call our model twisted quantum double 
model for reasons to be explained in Section TlVBI 



C. Equivalent Models 

Now that a 3-cocycle defines a twisted quantum dou- 
ble model, one may wonder since since a 3-cocycle rep- 
resents a whole equivalence class, whether two equiv- 
alent 3-cocycles, i.e., two representatives of the same 
equivalent class, define the same model. Let us con- 
sider two Hamiltonians Ho, a and HQ_ a i, respectively de- 
fined by two equivalent 3-cocycle a and a' that are re- 
lated by the 3-coboundary 8(3 of a normalized 2-cochain 
P : G 2 i-> U(l) that satisfy (3(x,e) = 1 = fi(e,x) for all 
x G G, 



a' (go, 91,92) 



S0(9o, 9i, 92)oc(jgo, 9i, 92 

P(9i,g2)P(go,gi92) 



P(gogi,g2)/3(go,gi) 



"(50,51,52), 



(13) 



where ^ £ G, and 8 is the coboundary operator. As 
each 3-cocycle is defined on three triangles (or equally a 
tetrahedron) such as in Fig [51 each 2-cochain /3 can be 
thought as defined on a triangle. Hence, Eq. IflSjl can be 
viewed as a local "gauge" transformation on a. 

We now check the relation between Hc, a ' and Ha, a - 
It suffices to check only the vertex operators A^(a') and 
A 9 v (a) because the face operators Bf have merely 8- 
functions as its matrix elements and are thus inert under 
the transformation in Eq. (|13p . Without loss of general- 
ity, we consider again the vertex operator on a trivalent 
vertex, as that in Eq. ©. By Eq. (flUj) , We immediately 
obtain the following. 



ff([23'],[3'4]) 



/3([12],[23'])/3([13'],[3'4]) 




(14) 



where the <5-function 5[3'3], fl is omitted for simplicity. The 
second term consisting of three a's is precisely the matrix 
element of A|(a). If we move the first fraction of (3 in 
the second equality of the above equation to the LHS, 
we readily see that the action of A^(a') on the rescalcd 
state 



/3([23],[34]) 




/3([12],[23])/3([13],[34]) 



matches perfectly the action of A^(a) on the original 
state. The above rescaling is clearly a local U(l) phase, 
which can be boiled down to the following local U(l) 
transformation on the basis states of triangles: 




i-> (3([ab],[bc]) e( - aM 




(15) 



where e(a, b, c) is a sign, which equals +1 if the enumer- 
ations a < b < c are clockwise on the triangle and — 1 
otherwise. In this new basis, A 9 (a 1 ) has the same matrix 
elements and thus the same spectrum as those of A 9 (a) 
in the old basis. 

There is a continuous deformation between any two 
3-cocycles related by a' = a.8/3. Define a 2-cochain 
P®(x,y) = P(x,y) 1 , with < t < 1, then a (t ) = aSfiW 
is equivalent to a for all < t < 1, with a' ' = a and 
a« =a'. The corresponding transformation in Eq. (|15p 
with (3 replaced by /3^ is a continuous local U(l) trans- 
formation; hence, there is no phase transition in the one- 
parameter family of systems with the the Hamiltonian 
H G a ( t ) from < t < 1. Thus we can conclude that the 
Hamiltonians Hq >cx i and Hq >01 due to two equivalent 3- 
cocycles a' and a indeed describe the same topological 
phase. 



III. TOPOLOGICAL OBSERVABLES AND 
SYMMETRIES 




AUa>) 



/([12],[23'I,[3'3])a'([23'],[3'3],[34]) 
a'([13'],[3'3],[34]) 




_ff([12],[23])/3([13],[34]) v a([12],[23'],[3 , 3])a([23 , ],[3 , 3],[34]) 



/3(123].[34]) 



X (L13'J,L3'3J,[34j) 



In Hydrodynamics, topological properties of fluid can 
be systematically studied by the diffcomorphism symme- 
try group acting on the fluid 2 ^. Similarly, the topological 
properties of a discrete model of topological phases can 
be systematically studied by the discrete version of diffeo- 
morphisms, which we shall call the mutation symmetry 
transformations of the graph. 

The symmetry we will be dealing with in this model 
are the mutations of the graph that preserve the spatial 
topology but not necessarily the local graph structure. 
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A Hcrmitian operator is a topological observable if it is 
invariant under the these mutation transformations. 

In most physical systems, the mutation (or diffeo- 
morphism) symmetry does not exist. Nevertheless, in 
the discrete models of topological phases proposed by 
Kitae v 30 ' 31 , and Levin and Wer>ii, the mutation transfor- 
mations to be constructed explicitly do have the space of 
the ground states of these models as invariant subspaces. 
Hence, we can use any topological observable, which is 
invariant under these mutation transformations, to char- 
acterize, at least partially, the topological phases in these 
models. One such topological observable is ground state 
degeneracy (GSD). 

In this section, we construct the mutation transforma- 
tions in our model and show that they are unitary sym- 
metry transformations on the ground states. Then we 
shall define and see, as an immediate consequence, that 
the GSD of our model is indeed a topological observable. 

All Bf and A v are mutually commuting projection op- 
erators, as proven in Appendix A. Thus the ground states 
are the simultaneous +1 eigenvectors of all Bf, A v . De- 
fine the ground state projection operator 



(l[B f )(l[A v ) 



fer 



and then the subspace of the ground states is 
H°={m\Pr\<S>) = \<S>)} 



(16) 



(17) 



Usually, symmetry transformations in a lattice model 
do not affect the lattice structure and are thus well- 
defined on a fixed lattice. The mutation moves in our 
model, however, take one graph to another. Since each 
graph r is endowed with a Hilbert space Tir and the 
Hamiltonian defined in Eq. (O, the Hilbert space and 
the Hamiltonian may be subject to changes under the 
mutation moves. 

It is known that we can always transform a triangular 
graph r to another one T' that triangulates the same 
Riemannian surface by a composition of the following 
elementary Pachner move o 32 ' 33 : 



Ti ■ 







(18) 



(19) 



(20) 



which are the generators of all mutation transformations. 
Each mutation generator fa : T — > V induces a linear 




transformation Tj : — > Hi 



T 



: ^2 "(fai^2], [V2V3), [V3V4}) 

[viV3]£LG 



To 



= ^2 a([qv 1 },[v 1 v 2 },[v2V 3 }) 

[v 1 q],[v 2 q], 

[v 3 q]&G 



T 3 



{.[viv 2 ], [V2V3], [V3V4]) 



(21) 




(22) 






(23) 



We now explain how we determine the linear properties 
of these operators. 

For Ti, we enumerate the four vertices in Eq. (|2ip by 
Vi < V2 < V3 < U4. The action of T\ does not change 
the degrees of freedom on the external edges, namely 
[^1^2], [1*21*3], [1*31*4]) and [U4U1], but only changes the 
[1*21*4] on the internal edge to [U1U3] on the internal edge 
in the new graph. The new group element [^11*3] runs 
over all group elements in G. The operator T\ also yields 

a V{\) phase a ([^11^], [U2W3], [W3U4]) I ' V2 ' 1 ' 3 011 tnc 
basis vector in the new Hilbert space. The exponent e is 
a sign function that assigns an exponent +1 or —1 to the 
3-cocycle a according to following rule. One first picks 
up either of the two triangles before the action of T\ , then 
notes down as a list its three vertices counter-clockwise, 
e.g., either (1*1, i*2> 1*4) or (vi,V3,v^) as in Eq. f2"Tj) . One 
then appends the remaining vertex to the list from the 
left, such as either (1)3, V\, 1)2, 1*4) or (vi, 1*2, 1*3, V4). If it 
takes even permutations to shuffle the list to completely 
ascending order, e = 1, which is the case in Eq. (I2ip . 
and £ = — 1 otherwise. 

As to T2 defined in Eq. (|22| . we suppose the order 
of vertices is v\ < V2 < 1*3. The action of T2 creates 
three triangles separated by three new edges that carry 
respectively three new group elements. We enumerate 
this new vertex by q, which is set to be less than v\ , such 
that the three new group elements are [qvi], [qv2], and 
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[qva], which are then averaged out in order not to enlarge 
the Hilbcrt space. 

The remaining factor in T2 is also a phase, which is 
in the form a ([uiWa], [^2^3], [vsq]) s ^ 9 ' Vl ' V2 ' V3 \ where the 
exponent is a sign depending on the orientation of the 
three triangles on the RHS of the equation. We determine 
the sign by first noting down the list of the three vertices 
clockwise from any of the three triangles of the basis 
graph after the action of T2, such as (q, 113,1)2) from the 
RHS of Eq. ([221) > then appending the remaining vertex to 
the beginning of the list, such as (i>i, q, V3, v%); if the list 
can be turned into ascending order by even permutations, 
e = 1, which is the case in Eq. ([22]) j an d otherwise 
e = — 1. In general, the enumeration q of the new vertex 
in Eq. (|22p can have any order relative the enumerations 
of the three old vertices; however, we assume q is the 
smallest therein for simplicity. 

As opposed to T2, T3 shrinks three triangles to a one, 
as in Eq. ([2"3"]l . the 3 -cocycle on the RHS of which is 

in general a ([V1V2], [^2^3], [v3Vi]) s ^ Vl ' V2,V3 , where the 
exponent is a sign depending on the orientation of the 
three triangles on the LHS of the equation. This sign is 
determined this way: One first reads off the list of the 
three vertices counter clockwise from any of the three 
triangles of the original basis graph, such as (l>2; 1*3,1)4) 
from the LHS of Eq. (|23j) . then appends the remaining 
vertex to the beginning of the list, such as (v\, v 2 , 1)3, U4); 
if the list can be turned into ascending order by even 
permutations, e = 1, otherwise e = — 1. To make life 
easier, in Eq. (|23|) . we consider only one case. 

Now we show that the mutation transformations gener- 
ated by T\, T%, and T3 are unitary symmetry transforma- 
tions on the ground states. In particular, T± is a unitarity 
of the entire Bf = 1 subspace of the Hubert space, in the 
sense that "Hp /_1 = Ti(Hp /_1 ). We denote the subspace 
of ground states of the Hilbert space Hr on a graph T by 
Hp- The proof consists of the following two steps. 

(i) . Mutation transformation preserve the space of 
ground states. 

That is, if T is a mutation transformation between 
two Hilbert spaces Hr to Hr', and if |$) 6 H r , then 
T|$) £ H T ,. 

It suffices to show that X!;-Pp = Pp/Tj for each mutation 
generator Tj and each state in H Bf=1 , where Pp and 
Pp, are the projectors respectively onto H r and H T , (see 
Appendix B). Note that however, we have H'T ^ TH in 
general. 

(ii) . Mutations are unitary on ground states. 

By unitary we mean: If T is a mutation transformation 
between two Hilbert spaces Hr to Hr', and if |$), \"^>) £ 
H r , then 

(T$ |T*) = ($ I*) . (24) 

It is sufficient to check Eq. ([23) for T X ,T 2 , and T 3 only, 
as seen in Appendix B. 

Consequently, there is a bijection between the ground 
states on any two graphs related by the mutation moves. 



Since two such graphs have the same spatial topology, 
the dimension of the ground state Hilbert space, i.e., the 
GSD of our model, is a a topological invariant and well- 
defined topological observable. Hence, our GSD can be 
taken as the trace of the ground state projector in Eq. 
(|16p . which is as we have just seen a Hermitian operator 
that is invariant under the mutations, namely, 

GSD = tr(P^), (25) 

where trace can be taken on any one of the graphs that 
are connected by the mutation moves, but the result is 
obviously independent of this choice. 

IV. THE GROUND STATE DEGENERACY 
AND TOPOLOGICAL DEGREES OF FREEDOM 

Ground state degeneracy (GSD) partially characterizes 
a topological phase. The nontrivial feature is that the 
GSD depends only on the spatial topology of the system. 
Two topological phases having different GSDs must be 
considered different. 

Another important characteristic of topological phases 
is the emergent fractional quantum numbers of the el- 
ementary excitations in these phases and the fractional 
statistics of the quasiparticles of these elementary exci- 
tations. The relation with GSD is that the GSD is equal 
to the number of species of the quasiparticles of the ele- 
mentary excitations. 

The significance of the GSD lies in the degrees of free- 
dom that are capable of distinguishing the degenerate 
ground states. The topological dependence of the GSD 
originates in that these degenerate degrees of freedom are 
global. An interesting question then arises: How to char- 
acterize these global degrees of freedom? Answering this 
question will enable us to (1) discern between two dif- 
ferent topological phases that have the same topological 
dependence of GSD and (2) understand better the rela- 
tionship between the global degrees of freedom in the de- 
generate ground states and the emergent fractional quan- 
tum numbers of the elementary excitations. 

In what follows we calculate the GSD of our model on 
a torus and then analyze the global degrees of freedom 
in the degenerate ground states. 

A. Ground state degeneracy on a torus 

The topological invariancc of the GSD of our model 
enables us to compute the GSD on the simplest triangle 
graph that triangulates the surface on which the model 
is defined. 

In the case of finite groups, the GSD of our model on a 
2-sphere is always unity because a 2-sphere has a trivial 
topology, in the sense that its fundamental group is triv- 
ial. This fact can be checked by following the approach 
to be presented shortly in this section. This is a common 
feature of all known models of topological phases. 
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A torus is the simplest closed surface with a non-trivial 
topology. Fig. 2] depicts the simplest triangle graph that 
triangulates a torus. 



2 9 
fa j ' 
1 



4 
/; 

9 3 



FIG. 4: Triangulation of a torus, with g, h 6 G. 

This graph has two triangle faces and only one vertex. 
But for the sake of assigning the 3-cocycles in A x easily, 
we use 1,2,3 and 4 to enumerate the sole vertex. This is 
perfectly fine because the boundary condition automati- 
cally merge the differently labeled vertices into one. We 
identify the boundary edge [12] with [34], and [13] with 
[24]. It is tricky to notice that the four enumerations 
can not be arbitrary. In Fig. U the orientations of the 
two boundary edges are consistently taken from higher 
enumerations to lower enumerations. 

The subspace H Bf=1 is spanned by the basis vectors 



{\g,h) \g,h eG,gh = hg} 



(26) 



corresponding to the assignment of [13] = g, [12] = h and 
[14] = gh = hg in the above graph. 

Since there is only one vertex in Fig. 01 we simplify 
the notation of A x at this mere vertex by A x , the action 
of which, according to its definition in Eq. Q , is 

A x | [13], [12]) 
= (a([13],[34'],[4'4])a([12],[24'],[4'4])-^ 

a({13'},{3'3],[34}y 1 a({12'},{2'2},{24'}) 
(a([l'l],[13'],[3'4'])a([l'l],[12'],[2'4'])-^ 
|[1'3'],[1'2']), (27) 

where the rule [l'l] = [2'2] = [3'3] = [4'4] = x holds 
for the new enumerations 1', 2', 3' and 4'. The coefficient 
consisting of the six 3-cocycles are determined as follows. 
We obtain (a ([13], [34'], [4'4]) a ([12], [24'], [4'4]) -1 ) at 

enumeration 4, then replace 4 by 4' to determine the 
rest of the factors. Next we obtain a ([13'], [3'3], [34])" 1 
at enumeration 3 and then again replace 3 by 3'. Hav- 
ing repeated similar steps at enumerations 2 and 1, we 
arrive at the above formula. One may derive a seem- 
ingly different coefficient by following a different path, 
e.g., 1 — > 2 — > 3 — > 4. But because of the topological in- 
variancc, the new coefficient can be brought to precisely 
the same as that in Eq. (|2"7| by applying 3-cocycle con- 
ditions, as one can check. 

Now we write down the action explicitly in terms of 



the group elements g, h and x. 

A x \g,h) 
=a(g, hx~ x , x)a(h, gx^ 1 , x)^ 1 
a(gx~ l , x, hx~ 1 )~ 1 a(hx~ 1 , x 



(28) 



a(x,gx ,xhx 1 )a(x,hx 1 ,xgx ) 1 \xgx ,xhx x ) 

One can verify the multiplication law A x A y = A xy by 
3-cocycle conditions ([3]) , which agrees with the result in 
Appendix El The ground-state projector is thus 



\G\ 



(29) 



Taking a trace of the ground-state projector (|2"9"]) com- 
putes the GSD, 

GSD =tr(|^£> s ) 



= ^2tgh,hg (g,h\A x \g,h) 

g,h 

=T77\ X! ^ghJig^hx^xg.gx^g^X' 1 , x) 
I ' h,g,x 

a(h, gx -1 , x)~ 1 a(gx~ 1 , x, hx~ 1 )~ 1 a(hx~ 1 , x, gx~ ) 
a(x, gx~ x , h)a(x, hx^ 1 , g)^ 1 (30) 

where the trace is evaluated in the subspace 7i Bf=1 . 

The seemingly complicated summation of the six 3- 
cocycles in Eq. (|30|) can actually be simplified in many 
cases due to a hidden simple mathematical structure. To 
see this, we shall first explore in the next subsection the 
algebraic structure in Eq. (|28j) . after which we come back 
to the simplification of the GSD. 



B. Topological degrees of freedom 

We now proceed to extract the algebraic structure in 
Eq. ([2"5]l and explore the classification of the topological 
degrees of freedom in the ground states, so as to reveal 
the deep mathematical significance of the GSD yet not 
fully discussed in the previous subsection. 

To this end, we rewrite Eq. (|28p as follows by applying 
appropriate 3-cocycle conditions (see Appendix C for the 
derivation) . 

A°\g,h) 

a(g, x" 1 , xhx~ 1 )a(x~ 1 , xhx~ x , xgx~ x ) 



a(x 1 ,xgx 1 ,xhx 1 ) 



a(h,g,x : ) 



\xgx 



\xhx~ 1 ). (31) 



a(g,h,x 1 )a(h,x 1 ,xgx *) 
By defining a new function in terms of 3-cocycles as 
def a(a, b, c)a(b, c, c~ 1 b~ 1 abc) 



Pa(b,c) 



a(b,b 1 ab,c) 



(32) 
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Va, b, c € G, and plugging it into Eq. (1311) . wc obtain 

P g (x~ 1 , xhx^ 1 ) 



A x \g,h) 



Pg{h,X-l) 



\xgx 1 ,xhx 1 ) 



where we define 

V g (h,x) 



ri 9 (h,x)\xgx ,xhx 1 ) 
(3 g (x~ 1 , xhx" 1 ) 



g (h,xrl) 



(33) 



(34) 



for any given g € G and h G Z g = {x 6 G\xg = gx}, the 
centralizer subgroup for g G G. Let h = g in the above 
definition, we have 



rj g (g,x) 



(3 g {x 1 ,xgx 1 ) 



(35) 



for all g, x € G, which can be quickly checked by di- 
rectly using the 3-cocycle condition. Interestingly, if 
x G Z g> h = {x G G|irg = gx,xh = hx}, we see that 
the U(l) number 



p 3 (h,x) =r) 3 (h, x) 



(3 a (x-\h) 

Pg^X- 1 ) 



(36) 



is actually a 1 -dimensional representation of the sub- 
group Zg : h C G. This is because p 9 (h,x)p 9 (h,y) = 
p 9 (h, xy), which is a consequence of A x A y = A xy on the 
ground states. 

It follows from Eq. (|33|) that the ground states are 
spanned by the vectors 



I ]^| zZ rfiKx) \xgx 1 , 



chx 1 ) 



g G G, he Z g 



(37) 



This tempts one to think that counting the GSD amounts 
to counting the elements in Horn (7ri(T 2 ), G) /con j, 
where the conj in the quotient is the conjugacy equiv- 
alence: (g,h) ~ (xgx -1 , xhx~ 1 ) for any x. This is in 
general not true, however, as one may over-count the 
states because in Eq. (|37|) . the terms that are summed 
over for some g and h may actually vanish, causing the 
corresponding states non-existing, as we now classify by 
studying the algebraic structure of the function /3 a de- 
fined in Eq. ([3"2]). 

Using the 3-cocycle condition of a, one can show that 
the function (3 a is in fact a normalized, twisted 2-cocycle 
satisfying twisted 2-cocycle condition, 



5(3 a (x,y, z) 



fix 1 ax 

(y,z)(3 a (x,yz) 



Pa(xy,z)p a (x,y) 
and the normalization condition 

/3 a (x,e) = f} a {e,x) = 1, 



(38) 



(39) 



for all a,x,y,z G G. The 8f3 a is called the twisted 3- 
coboundary of j3 a - 



Furthermore, when its variables are restricted to the 
centralizer Z a of a G G, (3 a clearly reduces to a nor- 
malized, usual 2-cocycle over Z a , which obeys the usual 
2-cocycle condition, 

Pa(y, z)Pa(xy, z)- 1 ^, yz)f3 a (x, y)- 1 = 1, (40) 

for all x,y, z G Z a . 

The function /3 a is closely related to the projective 
representations of Z a . In fact, each f3 a classifies a class 
of projective representations called ^-representations 
p : Z a — >GL {Z a ) obeying 



p{x)p{y) = /3 a {x,y)p(xy) 



(41) 



It is evident that the normalization condition corresponds 
to p(e)p(x) = p(x)p(e) = p(x), while the 2-cocycle 
condition in Eq. (f3"8")l corresponds to the associativ- 
ity p{x) (p(y)p(z)) = (p(x)p(y))p{z). In particular, if 
the 3-cocycles that define j3 a are the trivial one, i.e., 
oc = or = 1, then f3 a = 1, reducing p to the usual linear 
representations of Z a . 

In this paper, we are interested in the classification 
of /3 a ~representations of Z a with fixed a G G. Here we 
record a few important properties of this kind of repre- 
sentations. 

An element g G Z a is /? a -regular if fi a (g, h) = (3 a (h, g) 
for all h G Z a>g C Z a . Moreover, g is /3 a -regular if and 
only if all its conjugates are so, which can be verified by 
the 3-cocycle condition in Eq. ([3]). That is, g is /3 a - 
regular ^ [g] is /3 a -regular. We call [g] a f3 a — regular 
conjugacy class. In particular, each g is always @ g — 
regular. 

Let us denote all the conjugacy classes of G by C A and 
number of such classes by r(G). Since for any a, b G C A , 
Z a — Zb, it is convenient to denote these isomorphic cen- 
tralizers by Z A 7 obtained by any representative of the 
class. We henceforth collect any chosen set of repre- 
sentatives of all C A by simply R c = {g A G C A \A = 
l...r(G)}. 

For a G Z A , let the number of /3 a -regular conjugacy 
classes in Z A be r(Z A ,/3 a ). Clearly, we have 

r(Z A ,p a )<r(Z A ). (42) 

It is known that the number of incquivalent irreducible 
/3 a -representations of Z A is equal to r(Z A ,(3 a )- In par- 
ticular, in the case where /3 a = 1 because of the trivial 
3-cocycle a , we arrive at the familiar result that the 
number of the irreducible linear representations equals 
the number of conjugacy classes. Eq. (|4"2"]) states that 
irreducible /3 a -representations of Z a are fewer than the 
irreducible liner representations. 

The topological degrees of freedom is related to this 
classification of projective representations of Z A . To 
show this, we reexpress the GSD in Eq. (f3"0"f as 

GSD =— 5 gh , g h5gx,xg5hx,xhif{hi x ) 



g,h,x£G 



\G\ 

1 1 geGheZg xGZ S:h 



(43) 
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which can be further simplified, by the identity 



1 



E p 9 ( h ^) = 



1, his f3 g regular 
0, otherwise 



(44) 



where \Z g ^\ is the order of the subgroup Z g ^ with fixed 
g, heG. 

Here is the proof of Eq. (jlH). As shown below Eq. 
the phase p 9 (h 7 x) in Eq. ([M]) is a 1-dimensional 
representation of Z g ,h,\ it is the trivial representation 
p = 1 if h is /3 ff -regular and is otherwise a non-trivial ir- 
reducible representation (i.e., different from the identity 
representation). By the orthonormal condition 



1— E P s{j Hh,x) = S j , , 

J g.h\ ry 



(45) 



where j = corresponds to the trivial representation and 
j =/= a non-trivial irreducible representation, we obtain 
Eq. (0). 

Equation (gl} renders Eq. (J33]) as 



GSD = 



geGhGZ g 



\±aM y 
\G\ 

\Z A \ \Z n , 



g A M 



h€Z A ,A 111 1 



1, h is P g regular 
0, otherwise 



1, h is f3 g A -regular 
0, otherwise 



E 



r(Z 



A 



0- 



(46) 



In the last equality use is made of that |G|/|Z /l | = \C A \. 

According to the relationship between the number of 
/? s -regular conjugacy classes of Z g and the number of 
/^-representations of Z g as discussed above, the GSD 
can take the form 

GSD = E ^/^-representations of Z A ), (47) 

A 

where # stands for "the number of". 



C. Ground States Basis 

As promised in the previous subsection, we have sim- 
plified GSD evaluation in Eq. (|30|) to counting the rel- 
evant projective representations. Computing the GSD 
of our model on a torus amounts to counting the irre- 
ducible projective j3 g a -representations of each conjugacy 
class C A , then sum it over C A in G. 

Hence, the ground states on a torus can be labeled by 
pairs (g A , h) with g A running over Rq and h running 
over a set of j3 g a -regular conjugacy class representatives 
of Z A . Equivalently, the ground states can also be labeled 
by pairs (A, p) with A = 1 . . . r(G) and p labeling p^ , 
which are the irreducible (3 g A -representations of Z A . We 
posit that the basis vectors \A,p) can be defined as: 

1 



I A,p) 



E x%{h)\g,h), 



(48) 



where xfiW = tr/5^(/i) is the projective character defined 
as usual by the trace of the representation, and \G\ the 
order of G. Since the centralizers Z g arc isomorphic for 
all g in a conjugacy class C A , so are the set of irreducible 
^-representations of Z g for all g 6 C A . Therefore the 
same label p works for all ^-representations. We de- 
tail the construction of the isomorphism among the ir- 
reducible ^-representations for g G C A in Appendix [El 
The projective characters Xy.0 1 ) satisfy the following re- 
lation under simultaneous conjugation of g and h: 



XT (xhx-')=r 1 3(h,x)xf l (h). 



(49) 



for all x € G. Keep in mind that if h € Z A but 
h 7^ g, then h £ C A . Practically, for each conjugacy 
class C A with its representative element g A , if we find a 
/3 g a -representation p^ of Z A , we can construct the (3 g — 
representations for all other elements g of C A . Through- 
out this paper, we take the representations such that the 
relation (j4"9")) are always satisfied. 

Note that in general, the projective characters are 
not functions of conjugacy classes because the fact 
that ^(a)p£(&) = Pg{ a : fyp^iab) yields the relation 
Xlixhx- 1 ) = [P g {hx-\x)/P g {x,hx- 1 )]x%{h). Never- 
theless, the orthogonality and completeness relations of 
these projective characters still hold, namely, for all (3 g ~ 
regular elements a, 6 in G, 



7^7 E xlWxUh) = <W, 



hez„ 



\G J 



1, a conjugate to b 
"wise 



(50) 



where \Z g \ is the order of the subgroup Z g , and \C A \ is 
the cardinality of the conjugacy class C A containing a in 
the subgroup Z g . By Eq. ([SO)) , one can verify that the 
basis in Eq. (|4"8"|) is orthonormal. Moreover, if h is not 
/3 g -regular, then xfi(h) = 0, which is the very Proposition 
Q] proven in Appendix [E] 

We can now justify that \A, p) is indeed a ground state 
by its invariance under the action of the ground state 
projection operator P° defined in Eq. 

P°\A,I>) 

1 1 xEG 

= -7=£ E xlWih^xgx- 1 ^- 1 ) 
VI I xeG geC A 
hez g 

= E E XF^ixhx-^xgx-^xhx- 1 ) 



l G l 3 xeG„ eC 't 



hez„ 



geC A ,hez g 



VIM g > e c A ,h>ez g 



xeG 
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= -= X 9 fi ( h )\9,h) = \A,fx), (51) 

V l G l gec*hez t 

where Eqs. and (j4Uj) are used respectively in the sec- 
ond and third equalities, while substitutions g' = xgx~ l 
and h = xhx~ x are made to get the fourth equality but 
renamed back to g and h in the end. Therefore, we con- 
clude that the set of \A, fx) does furnish an orthonormal 
basis of the ground states, i.e., 

TL° = sp&n.{\A,n) : A=l...r(G), n = 1 . . . r(Z A , /3 gA )} 

(52) 

This uncovers the mathematical structure that clas- 
sifies the topological degrees of freedom in the ground 
states via representation theory. We start with our model 
specified by a 3-cocycle a over G and end up with the 
result that the topological degrees of freedom are deter- 
mined by the 2-cocycles (3 g over Z g . 

the ground-state basis vectors (A, //) label the set of 
all inequivalent irreducible representation spaces of the 
twisted quantum double D a (G), which plays a central 
role in the orbifolds by a symmetry group G of a holo- 
morphic conformal field theory We may understand the 
term "twisted" as twisting linear representations to pro- 
jective representations. We are not going to explain the 
details of the twisted quantum double, which is beyond 
the concern of this paper. But for completeness, we note 
here the multiplication law in the twisted quantum dou- 
ble D a {G): 

(P a ® x)(P b <8> y) = Pa(x,y)5 atXbx -i(P a ®xy), (53) 

for all a, b,x,y € G, where P a projects out a while x 
obeys the usual group multiplication with a projective 
phase factor. 

In particular, as to be shown in Section [Villi the un- 
twisted version of our model (i.e., when a = a ) turns 
out to be Kitaev's quantum double model (or, the toric 
code model), the GSD of which agrees with the num- 
ber of irreducible representations of the quantum double 
D(G) of the finite group G. Therefore, our model can be 
viewed as a deformation of the quantum double model 
by a twisting with the (3 a in Eq. (|32|) . which twists the 
linear representations of a group to the projective rep- 
resentations. This is mainly why we christen our model 
twisted quantum double (TQD) model. 



is known that the emergent fractional topological num- 
bers in the elementary excitations can differentiate such 
distinct topological phases. 

In this section, we shall first construct on the subspacc 
H B f =1 the topological observables then solve their eigen- 
problems to acquire the expected fractional topological 
numbers. These fractional topological numbers are re- 
lated to the fractional statistics of quasiparticles in the 
elementary excitations. Actually, there is believed to ex- 
ist a correspondence between the topological degrees of 
freedom in the ground states of the system on a torus and 
the local degrees of freedom of the quasiparticles in the 
elementary excitations. We shall come back to address 
this correspondence in Section HXl 



A. Topological observables as SL(2, Z) generators 

Consider the graph T on which the model is defined. 
In Section [TTTl we constructed the mutation transforma- 
tions that can change the local structure of the graph 
but preserve the graph topology, i.e., the topology that 
r triangulates. Under such mutations, the topological 
degrees of freedom of the ground states are intact. All 
such transformations are local. The ground-state projec- 
tor fl v A v can also be constructed from such mutations. 

Here, on the other hand, we look into the large trans- 
formations that alter the graph structure globally but 
still preserve the graph topology and lead to richer topo- 
logical observables. 

Again, since we are not interested in the local trans- 
formations of the graph, we need only to work on the 
simplest triangulation of torus as in Fig. 0J 

The transformations that change the topology are the 
familiar modular transformations, which form the group 
SL 2 {1i) that is generated by 

^(J-O 1 ). r-(Jl), (54) 

satisfying relations (ST) 3 = S 2 and 5 4 = 1. 

To cast the modular transformations in the form of 
3-cocycles, let us redraw the torus in Fig. [4] in the co- 
ordinate frame in Fig. [51 which illustrates the iS and T 
transformations on the torus. The S and T transforma- 



V. FRACTIONAL TOPOLOGICAL NUMBERS 

In the previous section, we studied the GSD as the 
simplest topological observable of our model. But topo- 
logical phases are only partially characterizes GSD. It is 
possible that two models specified by two inequivalent 3- 
cocycles have the same GSD but in the mean time, give 
rise to distinct topological phases. 

Hence, a natural question is how to differentiate two 
distinct topological phases if they bear the same GSD. It 



T: 



2^ I 

2 ^4 



FIG. 5: S and T transformations of a torus. 
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tions on the subspace are constructed as follows. 

We leave the details of the construction to Appendix 
ID] but claim here that 



=a([2'l] ) [13],[34])a([2'l],[12],[24])- 1 
xa([2'3'],[3'4'],[4'l])«([l'2'],[2'3'],[3'l])- 1 



xa([l'2'],[2'3'],[3 / 4' 1 



1' 



3' 



2' 



4' 



(55) 



where we set the order of the enumerations by 1' < 2' < 
3'<4'<1<2<3<4 such that the orientation of 
the two boundary edges are taken consistently. One sees 
that the wave function transforms oppositely in Fig. [5] 

Taking [12] = [34] = h, [13] = [24] = g, and [3'1] = 
[1'2] = [4'3] = [2'4] = x, S x is casted explicitly in terms 
of the group elements as 

S x \g,h) 

=a(xg~ 1 h~ 1 , g, ft,)a(xg _1 /i _1 , h, g)^ 1 
xa(xh~ 1 x~ 1 , xg~ x , g)~ 1 a(xgx~ 1 , xg~ x h~ x , h)^ 1 
xa(xg~ 1 h~ 1 x~ 1 7 xgx^ 1 , xg~ 1 )a(xgx~ 1 7 xg~ 1 h~ 1 x~ 1 , 
xa(xgx~ 4 , xg~ 1 h~ 1 x~ 1 , xgx~ 1 )\xh~ 1 x~ 1 , xgx^ 1 ), 

(56) 

where |[1'3'], [1'2']) d = Ix/i" 1 *" 1 , xgx~ x ). 
Similarly, we claim that T x behaves as 



T x 



=a([2'l],[13],[34])a([2'l],[12],[24])- 1 
xa([3"2'],[2'l],[13])a([4"2'],[2'l],[12])- 1 
xa([l"3"], [3"2'], [2'1]) a ([l"4"], [4"2'], [2'l]) _1 

4" V 



xa([l"3"],[3"4"],[4"2']) 




(57) 



where we set 1" < 3" < 4" < 2' < 1 < 2 < 3 < 4 as the 
order of enumerations, and when explicitly expressed in 



term of group elements, becomes 
T x \g,h) 

=a(xg~ 1 h~ 1 ,g, fyaixg^hT 1 , h, g)' 1 
xa(xhx~ 1 , xg~ x h~ x , g)a(xgx~ 1 , xg~ x h~ x , h)^ 1 
xa(xgx~ x , xhx^ 1 , xg~ 1 h~ 1 )a(xhx~ 1 , xgx~ x , xg^hr 1 )^ 1 
xa(xgx~ x , xhg~ 1 x~ 1 ,xgx~ 1 )\xgx~ 1 , xg~ 1 hx~ 1 ) (58) 

The s and T transformation arc defined by 

5 4E 5 ' T =^Y.r x - (59) 



\G\ 



The operators <S,T in Eq. (|55|) are a representation 
of the 5 and T matrices in Eq. (|5~i|) on the subspace 
T-L Bf=1 of the model. Indeed, a direct evaluation by the 
3-cocycle conditions verifies (ST) 3 = S 2 and S 4 = 1. 

On a torus, the vertex operator A x and the modular 
transformation operators S x and T x comprise an inter- 
esting algebraic structure, namely, 



s x A y = S xy = A x s y^ 

T x A v = T xy = A x T y 



(60) 



We shall not prove this here, as it can be done straight- 
forwardly by manipulating the 3-cocycles in the above 
equations. This algebraic structure results in the follow- 
ing important reexpression of the S and T operators. 

p°T=T = T e P°, P°S = S = <S e P°, (61) 

the first of which is proven as follows. 



t= — Yt x 

\G\ Z - 



Sq.JHJ 1 



=T e P° = P°T e , 



where e is the identity element of G. The proof of the 
second relation in Eq. (fBTjl follows likewise. This indi- 
cates that the operators S and T are indeed topological 
observables and symmetries in H°. 

We can lay the ground states in the basis composed of 
the eigenvectors {&k} of T, 



(62) 



where Ok is a t/(l) phase, and k = 1,2,...,GSD labels 
the degenerate ground states. These eigenvectors will be 
identified with \A, yt) in the next subsection. 

We remark that T also has other eigenvectors, whose 
eigenvalues are zero, which is implied by the first relation 
in Eq. (foTj) . These zero eigenvectors are actually the 
excited states of the model; however, we are not going to 
dwell on them in this paper. 

Hence, one can regard the eigenvalues 6}. of T as a set 
of topological numbers of the model. In fact, from Fig. [5J 
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T can be viewed as a global twisting of the system, and 
thus its eigenvalues 9k can be regarded as the topological 
spins of the topological sectors \<&k)- 

Another set of topological numbers are the 5-matrix 
of the topological sectors, 



(63) 



where i, j = 1, 2, GSD. This matrix is orthonormal: 



j 



= s. 



(64) 



Above all, apart from GSD, we obtain two more sets 
of topological numbers, {Ok}, and {sij}, to characterize 
the topological phases in our model. 

We remark that we have presented here a novel deriva- 
tion of the modular S and T matrices, which is purely 
based on our model and in terms of 3-cocycles of G, with- 
out resorting to any theory of group representations. 



B. S and T Matrices 

We now offer concrete solutions of the topological num- 
bers {9k, s^}, which are tied to the projective represen- 
tation theory. We emphasize that the topological ob- 
servables S and T are defined on the subspace T-L Bf=1 , 
whereas the solutions to their eigen-problems are to be 
obtained on U° C U B t =1 . 

In the following we diagonalizc the T matrix in Eq. 
()59p . One should bear in mind that the transformation 
T acts non-vanishingly only on the ground states. In 
Section IIVC1 we see that the ground states are spanned 
by the orthonormal basis \ A, fx) defined in Eq. (|4"5)l . with 
A running over all conjugacy classes of G and fx over the 
irreducible (3 g a -representations of Z A . It turns out that 
| A, li) arc eigenvectors of T as we demand, Here we sketch 
the proof. In the \A, p) basis, the action of T becomes 

T\A, fx) 
=T e P°\A,fi)=T e \A,fi) 

=W\ £ xim^g-^A,") 

1 1 sec" 



xl (g A ) 

dim,, 



\Ap), 



(65) 



where dim^ is the dimension of the representation fx, the 
second row uses Eq. (|6"Tj) and Eq. (|5"Tj) . and in the fourth 
equality use is made of the inverse transformation 



Therefore, the basis vectors \A, fx) are indeed the eigen- 
vectors of T, with the eigenvalues 



9 A 



xj {g A ) 

dim,. 



(67) 



Clearly, the projective characters are the ground-state 
wave functions of the system, collapsed in the basis vec- 
tors that are the eigenvectors of the T matrix. 

The above calculation guides us to interpret 9 A from 
the representation theory as an invariant that charac- 
terizes the representation ~pF . More precisely, for any 
g G C A , the matrix p^Ag) commutes with all other ma- 
trices p^(/i) for h G Z g , as 

WW = fif^WW = f%(hK(9), (68) 

where the second equality can be checked by a direct 
evaluation in terms of 3-cocycles. From Schur's lemma, 
the matrix p^(.g) is a multiple of the identity matrix 



dim u 



1. 



(69) 



Moreover, by setting h = g in Eq. (|4*9^1 and using Eq. 

(|3"5j) . we find that x X9X (^ff^ -1 ) = X 9 (ff) ; indicating 
that the topological number 9 A = x^(3)/dim M is indeed 
an invariant on C A , associated with the representation 
Pi of Z A . 

The topological spin 9 A are given in terms of 3- 
cocycles a by 



A\PA 



with 



UJA 



= f ' [ a(g,g n ,g) 



PA — 1 



(70) 



(71) 



for conjugacy class C A of G, and pa the degree, i.e., the 
least integer such that g PA = e, where g n is the power of 
g. The oja is independent of the choice of g G C A and 
thus a conjugacy class function. The relation is verified 
by applying Eq. (JUJ) to (p A ) PA = (9 A ) PA l. Therefore, 
the topological spin 9 A takes values in the Pa~^ roots 
of ui a- Moreover, each of the pa distinct p^-th roots 
appears precisely r(Z , (3 9 )/pa times in {9 A } for all fx, 
where the ratio r(Z A , f3 9 )/pa is an integer 44 . 



\g,h) = 



i 



r{Z B ,fl g ) 



E xl{h)\B,v) 



(66) 



which is defined in W° only, with g G C B is assumed. 
Appendix [El proves Eq. ([r?5|) step by step. 



To see this, one observes that each element in Z g can be uniquely 
written as g n h with n = 0, l,...,^ for some h G Z g , and 
there are pa 1-dimensional representation of Z g by ffl(g n K) = 
cxp(2irijn/pA)- Then for each representation fi there exist 
a fi' such that p^,(g n h) = pi (g n h)pfi (g n h) , and thus that 
S A , = exp(27rin/p A ) 
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Similarly, the 5-matrix can also be evaluated in terms 
of the projective characters. We record as follows the 
final formula for the 5-matrix while detail the proof in 
Appendix [E] The 5-matrix reads 



s (An)(Bv) = (A,(i\S\B, 



\G\ 



xUh)x h Aa)- (72) 



g£C A ,heC i: 
gh=hg 



Again, we take the projective representations ]fi such 
that the projective characters are related by Eq. (|4"9")l . 
This general result of the 5-matrix actually offers an 
answer to one of the open questions listed in Ref2£. 

The mathematical significance of the S and T- 
matrices is the following. They are the invariants carried 
by the projective representations of Z A for all conjugacy 
classes C A , in which the (3 g functions play a crucial role. 
All these ingredients are well-organized by the represen- 
tation theory of a twisted quantum group (or, a twisted 
Hopf algebra) , called the twisted quantum double D a {G) 
of the finite group G, which is parameterized by a 3- 
cocycle a. All the irreducible j3 g A -representations for all 
conjugacy class representatives g A form the linear irre- 
ducible representations of D a (G). 

Usually, the irreducible representations of a (twisted) 
quantum group classify the anyonic quasiparticle species. 
The invariants of each irreducible representation identi- 
fies the fractional topological quantum numbers of the 
corresponding quasiparticle. The 5-matrix has the origin 
as a braiding operation that exchanges any two of these 
quasiparticles, while the 7"-matrix contains the statis- 
tical spins of the corresponding quasiparticles which are 
determined by the braiding operation. For the discussion 
of the S and T-matrices for the twisted quantum double 
D a (G), see Ref^^. 

In this section, we have reproduced from our S and T 
operators in terms of 3-cocycles the familiar S and T— 
matrices in terms of projective characters for the twisted 
quantum double D a {G), which were originally obtained 
from representation theory, according to a braiding oper- 
ation. Our calculations are carried purely in the ground- 
state subspace and root in the large transformation of the 
spatial graph of the system on a torus. We expect that 
the quasiparticles in the elementary excitations will be 
classified by the same topological numbers in the way 
that the GSD equals the number of the quasiparticle 
species, the S and T-matrices on the ground states are 
the same as those of the quasiparticles. 



VI. TOPOLOGICAL NUMBERS AND 
TOPOLOGICAL PHASES 

We believe that the topological phases are classified by 
the topological numbers {GSD, 9 A S(Au) (Bv)}- In ai l ex- 



amples discussed in Section fVIfl we observe that they are 
classified by the third cohomology classes of a, i.e., any 



two models Ha, a an d Ha, a ' have the same topological 
numbers if and only if a and a' . 

In this section, we study how the topological numbers 
depends on the cohomology classes of a. 



A. When 3— Cocycle is cohomologically trivial 



In Section Hi C| we have shown that two equivalent 3- 
cocycles define equivalent twisted quantum double mod- 
els, which consequently should describe the same topo- 
logical phase. We now study this topological phase in 
more details. 

We begin with a special case, where the 3-cocycle of 
our model belongs to class of the trivial 3-cocycle a . 
Such a 3-cocycle can take the form of a 3-coboundary: 



oc(x,y,z) = 6P(x,y,z) 



fi(y,z)P(x,yz) 
f}(xy,z)f){x,y) 



(73) 



where f3{x, y) is any normalized 2-cochain, i.e., any func- 
tion GxG-> U{\) that satisfies /3(e,x) = 1 = (3(x,e) 
for all x £ G. To be seen in Section IVIII1 such a model 
is equivalent to Kitaev's quantum double model. 

The corresponding twisted 2-cocycle j3 g , defined in Eq. 
(|32jl , is automatically trivial such that it has the freedom 
to be written as a twisted 2-coboundary: 



where 



/3 g (x,y) = 6e g (x,y), 



e g (x) = fi(x,x 1 gx)/3(g,x) 1 



(74) 



(75) 



is a twisted 1-cochain, whose twisted 2-coboundary 
reads 



Se g (x,y) = e g (x)e g (xy) 1 e x -i gx (y), 



(76) 



for all g, x, y 6 G. From the relation ([75]) . we inevitably 
notice the following constant : 

e g {h)e h {g) = 1 and e g (g) = 1, Vg,heG. (77) 

This is indeed a constant as it is clearly independent of 
which a is picked in its equivalent class. 

By the form of j3 g in Eq. (f87f . the irreducible f3 g - 
representations of Z g are in one-to-one correspon- 
dence to the irreducible linear representations p M , by 



(78) 



for all he Z g , by which one can directly check the defi- 
nition property (|41|l . 

The ground states in Eq. (|4"8j) now become 



\TGl E z a WMg,h), 

V 1 1 g£C A .heZ„ 



I Am) 

where xf± = t r P^ is the usual character. 



(79) 
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Eq. (|T6"|) implies that /3 g (x,y) = P g (y,x) for all x,y & 
Z g with xy = yx. This means that all elements in Z g are 
/3 g -regular and that rj g (h,x) = 1 for all x £ Zg.h- Hence, 
the GSD in Eq. gSJ) now reads 



GSD = E E E 



1 



Hom(7n (T 2 ),G) 



con] 



(80) 



where the quotient means the equivalence (g, h) ~ 
(xga; -1 , xhx^ 1 ) for any x £ G. 

By ([75)1 and the constraint (f77)l of the e n , the topolog- 
ical numbers and the 5-matrix are expressed by 



9 A 



xl (g A ) 

dim,, 



xj (g A ) 

dim,, 



and 



\G\ 



E xl{h)x h M- 



gec A ,heC E 

gh=hg 



(81) 



(82) 



When the 3~cocycle is a = 1, the ground states are 
labeled by the usual irreducible linear representations of 
all the centralizers Z A C G. For a £ [a ] but a ^ a , the 
ground states are labeled by projective representations, 
which arc related to the corresponding linear represen- 
tations by merely a phase, of all the centralizers, since 
fig 7^ 1; however, all topological numbers are the same as 
those in the case of a , as they should be. 



B. When twisted 2— cocycle is cohomologically 
trivial 



When the 3-cocycle a ^ [a ], it could still be "trivial" 
at a lower level, in the mathematical sense that the 2- 
cocycle j3 a it defines in Eq. (f3"2"j) is cohomologically trivial, 
i.e., this p a is actually a twisted 2-coboundary: 



Pa(x,y) = 5e a (x,y), 



(83) 



for all a,x,y £ G. Note that however, the twisted 1- 
cochain e a in this case does not necessarily have the 
closed form in Eq. (|75|l in general because a is not co- 
homologically trivial; hence, e g (g) ^ 1 in general. The 
twisted 2-cocycle condition in Eq. f3"8"j) yields 



n°{h,x) 



e xgx -i(xhx 1 ) 



(84) 



for all h £ Z g and x £ G, which is unity for all x £ Z g ^. 

Similar to the previous case, the ground-state subspacc 
in the current case are also spanned by the basis vectors 
\A,(jl) of the form in Eq. (|75|) . where n labels the p g - 
representations of Z g , which are again related to the 
usual linear representations p M by Eq. (|78|) . 

Since Eq. ([M]) renders all elements in Z g /3 g -regular, 
as before, the GSD in this case copies that in Eq. (fBUl) . 



By the form (|75|) , the topological numbers 0£ and the 
s matrix are related to e g (h) by, 



9- 4 - f Jn A \ X ^ — 



and 



(85) 



(86) 



gh=hg 



The GSD is the same as the one in the [a°]-modcl, 
a result of the cohomologically trivial (3 g . Nevertheless, 
the topological numbers 8 A and the 5-matrix charac- 
terize the difference between the current model and the 
untwisted model. The Eq. ([55]) exhibits the physical rel- 
evance of the phases e: They endow each basis ground 
state \A, a) an extra spin factor e g A(g A ) in addition to 
the topological spins in Eq. (f8"l"T) . 



C. General properties of the topological numbers 

Since our model is defined in terms of 3-cocycles a, all 
of the topological numbers discovered must depend on a 
at the bottom level. Notwithstanding this, the special 
cases belabored above imply that in general, some topo- 
logical numbers may have a higher-level dependence of 
a via the quantities derived from a, such as the equiv- 
alence class of a and the twisted 2-cocycles. We now 
list as follows the general characteristic properties of the 
topological numbers {GSD, 8 A , S(An){Bv)} of our model 
on torus. 

1. The set {GSD, 6 A , s^^g^} of all topological 
numbers depends on the equivalence class [a] £ 
H 3 (G, U(l)) of the a that defines the model. 

2. The GSD depends only on the equivalence classes 
[p gA ] £ H 2 (Z A ,U(l)) for g A £ R c , A = 
1, . . . , r(G), independent of the representatives g A . 

3. The topological spins {9 A } are classified by 
{r(Z A , Pa), wa} for all conjugacy classes C A . 

We now elaborate on the two properties above in order. 

For property 1, one can check that any two equivalent 
3-cocycles a' and a related by a' = aSP, as in Eq. (fHij) . 
give rise to two equivalent twisted 2-cocycles P' a and p a 
related by a twisted 2-coboundary as follows. 



p' a (b, c) =0 a {b,c)6e a (b,c), 



(87) 



where the twisted 2-coboundary Se happen to be those 
defined in Eq. (|75|) and Eq. (|76| . The constraint in 
Eq. C2| implies that GSD' = GSD, 6' A = 6 A and 



S (Ap,){Bu) — s (A^){Bu) 



Property 2 is manifest in Eq. (|4"6"1) . in which the GSD 
is a sum of the numbers r(Z A , P 9 a) over all conjugacy 
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classes C A of G. One further confirms this by looking 
at Eq. (Jl3"|) . where the GSD is a sum of the phases p 9 
as a function of (3 g defined in Eq. which is a 1- 

dimensional representation of Z g> h with h € Z g . Suppose 
two 3-cocycles a and a', not necessarily inequivalent, 
that define (3 a and j3' a equivalent for all a £ G, in the sense 
that there exists functions e a ■ G — > U(l) parameterized 
by a £ G, which satisfies e a (e) = 1, such that f3 a and f3' a 
arc different by merely a twisted 2-coboundary, namely, 



e a {x)e x -i ax (y) 



Pa(x,y), 



for all a, x, y £ G. In the restriction to x £ Z a and 
y £ Z atX , one sees that and j3 a are equivalent up to a 
usual 2-coboundary over Z a . By Eq. ([55]) and Eq. (03]), 
we obtain 

p' 9 (h,x)=p 9 (h, X ), 

which verifies that GSD' = GSD. Furthermore, Eq. (gj) 
states that the sum of p 9 (h 7 x) over Z g> h is either unity 
or zero, regardless of which representative g £ C A is 
chosen. Therefore, property 1 holds, as expected from 
the analysis in Section ["II CI 

Property 3 is straightforward. As discussed in Sec- 
tion [V] each of the pa distinct pa -roots appears pre- 
cisely r(Z A ,f3 9 )/pa times in {8 A } for all p. Then 
uj a = uj a i yields 9 A = 8' A up to a relabeling of p = 
0, 1, r(Z A , p A ) - 1, assuming r{Z A , p A ) = r{Z A ,P' A ). 



The m cohomology classes of 3-cocycles are generated 

a : Z m x Z m x Z m t/(l) 

a(a, 6, c) = exp |^2 a [fo + c - (& + c))]| , 

where a, 6, c € Z m , and (a) is the residue of a mod m. 
By "generated" we mean that the m classes of 3-cocycles 
can be represented by the powers of a in Eq. (|89p as 



{a fe (a, 6, c)| fc = 0, 1, ...,m- 1}. 

One verifies that a m = 1 and a k a l = a k+l . 
The (3 g for each a fe has the form (f53"|) . with 



2ttI , 



e (&) = exp{ — x a x b} 



(90) 



(91) 



Hence each a k gives rise to GSD = to 2 . 

The linear characters are Xni a ) = exp(27ri/ia/TO,), for 
p = 0, 1, ...,to — 1. Applying this to Eqs. (|55|) and 
yields 



^ = exp [27ri(fca 2 + TOa/x)/m 2 ] 

s ( at i),(biy) = — exp {-2ni{2kab + m(av + bp)]/m 2 } 



(92) 



B. G = Z 2 and Z 3 



VII. EXAMPLES 

In this section, we explicitly compute various exam- 
ples of our model, making contact with the structure dis- 
cussed in the previous sections. We declare that some 
parts of the results here are adapted from certain known 
results ^ 28 ' 35 of 3-cocycles and projective representa- 
tions that were otherwise discovered in studies of con- 
formal field theory by means of representation theory, 
which now, however, as we show, become applicable to 
describing topological phases, owing to the lucid connec- 
tion revealed by our model between topological phases 
and group cohomology. In all examples to be discussed 
in this section, any two models Hc, a and Hc, a ' have the 
same GSD and satisfy the condition in Eq. ([77]) if and 
only if a and a' are equivalent via Eq. (TT3|) . 



A. G = Z m 

When G is the cyclic group Z m of order m, it is known 
that the cohomology group is H 3 (Z m , U(l)) = Z m , and 
hence there are to inequivalent classes of 3-cocycles 3 ^. 

We denote by a £ {0, 1, to — 1} the elements of Z m . 
The multiplication in Z m is a ■ b = a + b mod to. 



In the special case where G is the simplest finite group 
Z2, there are two classes of 3-cocycles. The first is the 
trivial one, namely, 

a°(a,6,c) = l, (93) 

where a, 6, c = 0, 1 are elements of Z2. 

The second one is given as follows, according to Eq. 

(83), 



a 1 (a, b, c) = 1, for all other a, b, c 

We recognize the a model as Kitaev's toric code 
model, or dual to the Levin- Wen model with the 6j 
symbols determined by irreducible representations of Z2 . 
The a 1 model is dual to the Levin-Wen model with 
the 6j symbols determined by the semisimple irreducible 
representations of the quantum group U q (sl(2, C)) for 
q — exp(i7r/3), up to a local unitary transformation (see 
Section [X| . The form the complete solutions to Levin- 
Wen models with Z2 fusion rule. The topological spins 
8^ for ground states (^) are given in Table fl] 

When G = Z 3 , there are three classes of 3-cocycles, 
denoted by a , a 1 , and a 2 . The a model is dual to 
the Levin-Wen model with the 6j symbols determined 
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(8) (?) (o) (i) 



1 i i 

1 1 i 



TABLE I: 9f, for models with G = Z 2 





(8) (?) 


(a) 


(o) (1) 
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i) © 


(f) (1) 
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i i 


1 


1 e* 




¥ i 


2tt1 j»i 

e 3 e 3 


a 1 


i i 


1 


2tt> 8xi 


e~ 


4iri SttI 

"5" e"S" 


27ri 4vri 

e~$~ e "5" 


a 2 


i i 


1 


47ri 8iri 

e s e s 


e~ 


27ri 27ri 
9 e 9 


87ri 4iri 

c 9 e 9 



TABLE II: 



for models with G = Z3 



has the form 



with 



(4 1) ) a ( 6 ) = ex P{^Jai x M, 

(4 2) ) a ( fo )-cxp{^a 2 x6 2 }, 

(4/ 2) )a( & ) = exp{^ ai x 6 2 }. 



(96) 



Therefore the associated /3 Q for all 3-cocycles are 
equivalent, and correspond to the trivial element in 
H 2 (1^,1/(1)) = Z m , though the second cohomology 
group itself is non-trivial. 

We conclude the models specified by all 3-cocycles 
have GSD = m 4 . 



D. 



G — X Zirn X 



by the irreducible representations of Z 3 (or equivalent to 
Kitaev's model with Z 3 ), whereas the a 1 and a 2 models 
are not dual to any Levin- Wen models with the Z 3 fu- 
sion rule, as there is only one Levin- Wen model in the 



circumstance. The topological spins 
lated in Table HH 



are also tabu- 



Consider the complex conjugation JC. The models 
H% 2 a 1 H^.a 1 an< i H% 3 a o are invariant under JC, whereas 
the models H% 3 , a i and Hz 3>a 2 ar e not. The complex con- 
jugation JC transforms the topological spins in the H^ 3 a i 
model to those in the -ffz 3 , a 2 model, as seen in Table HI1 
The flz 3jQ i and H% 3t0t i models are the simplest mod- 
els that break the complex conjugation symmetry which 
persists in the Levin- Wen models. 



G — Zm X Zirj 



The simplest non-cyclic Abelian group is G = Z^. 
But nothing is really new compared with the case where 
G = Z TO . The cohomology group # 3 (Z^, U(l)) = Zf n 
has three generators. We will label the group elements in 
G = Z^ as pairs a = (01, a 2 ) with a\, a 2 = 0, 1, m — 1. 
The multiplication is the obvious one (ai, a 2 )(&i, 6 2 ) = 
(ai + a 2 mod m, b± + 6 2 mod m). The three 3-cocycle 
generators are 

af\a,b,c) = exp{^ai(6i +c\ - (h + ci))}, 



f\a, b, c) = cxp{— a 2 (b 2 + c 2 - (62 + c 2 ))}, (95) 
2tti 

(a, 6, c) = exp{— ai(6 2 + c 2 - (6 2 + c 2 ))}, 



„( 12 )/ 



where (x) = x mod m is the residue of x. The m 3 classes 
of 3-cocycles are the products of powers of these three 
generators. The j3 a function for all these three generators 



When it comes to the case of G = 1? m , things become 
more interesting. We label the group elements by triples 
a = (ai, a 2 , a 3 ) with a l7 a 2 , a 3 = 0, 1, ...,m— 1. The coho- 
mology group i? 3 (Z 3 ra , U(l)) = l? m has seven generators, 

a { j\a, b, c) = exp{^ai(6j + Cj - (b t + Cj))}, 

aff (a, b, c) = exp{^o i (6 j + cj - (bj + Cj ))}, (97) 

auiia, b, c) = exp{— a!6 2 c 3 }, 
m 

where 1 < i < 3 and 1 < i < j < 3 are assumed respec- 
tively in the first two lines, and (x) is the residue of x 
mod m. 

The P a function for the first two types has the form of 
Eq. (SSI, with 



(e ( f\{b)=eM^a j xb j }, 

(4/ fe) )a( 6 ) = CX P{5 a J" X M> 



(98) 



But am cannot be decomposed as in Eq. (f8"3"j) . This 
provides further classification of the models defined by 
the m 7 3-cocycle representatives. Precisely speaking, 
each 3-cocycle class has a representative of the form of 
a product of powers of the seven generators. They are 
classified into m classes, depending on the power of am 
they contain, namely the form a q jjj for q = 0, 1, m— 1. 
The 2-cocycle for each a^jj is given by 

27ri 

f) a (b, c) = exp{ q(aib 2 c 3 - 6ia 2 c 3 + &ic 2 a 3 )}. (99) 

n 

If a\ and a 2 belong to two different classes with the pow- 
ers qi and <j 2 7^ q\ , the corresponding models have inequal 
GSD in general. 

Then the GSD in Eq. for the model specified by 
a depends on the class only, i.e., the power q of am. 
Specifically, the GSD is determined by Eq. P3"|) , with f3 a 
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9 = 


9 = 1 


q = 2 


g = 3 g = 4 


rn 


= 2 


64 


22 






m 


= 3 


729 


105 


105 




m 


= 4 


4096 


400 


1408 


400 


m 


= 5 


15625 


745 


745 


745 745 



TABLE III: GSD with G = Z m xZ ra xZ„ 



For n > 3, things arc similar to G = H? m case. The second 
and the third cohomology groups are 

H 2 (z^ l ,c/(i))^z^- 1 )/ 2 , 

ff 3 (Z™ , 17(1)) ~ Z"+"(™-l)/2+n(n-l)(™-2)/3!. ( 102 ) 

Like in the l? m case, there are three types of 3-cocycles 
taking the following form 



given by Eq. (|99[) . According to the analysis in Section 
IIV Al computing the GSD amounts to count the number 
of /3 a -regular conjugacy classes of Z a , for each conjugacy 
class representative a of G, then sum over all a. The 
GSD is 



GSD = j T J[[(p k * ~ m+p- 1 + P - 2 ) + 1] , (100) 



where / = m/gcd(g, m) is the greatest common divisor 
of q and m, p and k p are the prime number and the corre- 
sponding power in the prime decomposition / = Y[ p P kp ■ 
The GSD in m = 2, 3, 4, 5 cases arc given in Table IIII1 
which can be computed cither from the original formula 
Eq. fl-23) or the final reduced formula Eq. {IDD]) . 



27ri 

a { j'(a, b, c) = exp{ — + c. t - (b { + Cj))}, 



4 l) ( 



aW (a, 6, c) = exp{ — -ai(6j + 

?7Z 



2ttI 

"2 1 

( - 



(6j+ty))}, (103) 



-2?ri 



ayf, (a, 6, c) = cxp{ a^c*,}, 

in 



where 1 < « < 3, 1 < i < j < 3, and l<i<j<fc<3 
are assumed respectively in the three lines, the num- 
ber of three types of generators are n, n(n — l)/2, and 
n(n— l)(n— 2)/3!, corresponding to the number of gener- 
ators in H 3 (Z^ n , U(l)). Topological phases are classified 
by the elements in i7 3 (Z^, U(l)). The 2-cocyclcs (3 a ob- 
tained from type-Ill generators correspond to non-trivial 
elements in i7 2 (Z™ , U(l)). 



The models with q = 0, i.e. with only a/ and an as the 
defining 3-cocycles, possess similar topological numbers 
{GSD, 8, s} as those in the previous examples for G = Z m 
and G = Z m x Z m . The 9 and S are derived from the 
linear characters 



2?ri . 

Xn\. x ) CX P 1 — (Mi^i + ^2X2 + m x 3) 



(101) 



with /ii,^2,M3 = 0, l,...,m— 1 labeling the irreducible 
representations of G = Z m x Z m x Z m . 

The models with q =/= 0, i.e., involving a/// in the defin- 
ing 3-cocycle, possess more interesting topological num- 
bers. Though the finite group G is Abelian, the topolog- 
ical charge of the ground states are non- Abelian. 



F. G = D m 

The simplest non-Abelian groups are the dihedral 
group D m . Specifically, D3 (equivalent to the permu- 
tation group S3) is the simplest non-Abelian group. We 
will only consider odd m here, in which all 3-cocycles 
can be decomposed as in (|83p . 

We will label the elements in D m by pairs (A, a) for 
A = 0, 1 and a = 0, 1, m — 1. The multiplication law 
takes the form 

(A, a)(B, b) = ((A + B) 2 , ((-l) B a + b)J (104) 



E. 



G — 



Any Abelian group is isomorphic to ZJ^ for some inte- 
ger m and n. We have investigated the cases with n < 3. 



where (x) 2 = x mod 2 and (x 



= x mod m means tak- 
ing the residue. 

The cohomology group H 3 (D m , U(l)) = Z 2m has only 
one generator of the 3-cocycles: 



a((A,a),(B,b),(C,c)) = exp <j 



(-l) B + c a [(-ifb + c ((-Ifb + c) m ] + ^ABC 



(105) 



The representatives of each 3-cocycle class takes the The (3 a for this 3-cocycle generator takes the form (|83jl , 
form a p for p = 0, 1, 2m — 1. 
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with 



e ( A, o) ((#,&)) 

-xp(g{fe[(-l) B a + 2A6]- Ab* + ^AB 



(106) 



for all (A, a) and (B,b) that satisfy (A, )(£?,&) = 
(B,b)(A, a). The GSD is the same for all 3-cocycles, 
and are given by GSD = m + 7 , 



VIII. KITAEV'S QUANTUM DOUBLE MODEL: 
a IS TRIVIAL 



In this section, we show that in the special case where 
the 3-cocycle is trivial, our model becomes the Kitaev's 
quantum double(QD) model. By "trivial" we mean that 
the 3-cocycle takes the constant value 1, 



a°(x,y,z) = l, for all x, y, z G G 



(107) 



With this a , the definition © of A v operator is re- 
duced to 



A, 





(108) 



Then the model defined by © and (JT|) becomes the 
familiar Kitaev's quantum double model on triangle 
graphs. With a nontrivial 3-cocycle a, our model can 
be viewed as the twisted version of Kitaev's QD model, 
where the twisting is specified by the 3-cocycle a. We will 
explain the twisting in more detail in the next section. 

To gain more intuition, we would like to briefly review 
Kitaev's quantum double model in the language of gauge 
theory. 

To set up a gauge theory on the graph T, we need to 
specify the connections and the gauge transformations. 
Each basis vector in Eq. (QJ corresponds to a connection, 
namely, an assignment g : E — > G to each edge e of T 
a group clement g e of G. A gauge transformation h on 
F is an assignment to each vertex v a group element h v 
of G. The action L(h) of a gauge transformation fcona 
connection g is given by [L(h)g] e = 3t( e )3eff7(e) for cach e > 
where s(e) and t(e) are the starting and ending vertices 
of the edge e. For example, on one edge e orienting from 
v\ to V2, the action of a gauge transformation h is 



L(h): 



vi v 



„i j . v 2 ) (109) 



The action of any gauge transformation can be decom- 
posed into local operators defined at each vertex. We 



denote by L v (h v ) the action of a local gauge transforma- 
tion of at vertex v, which is defined as, 



L v (h v ) 




The Hamiltonian of Kitaev's QD model is 



/ 



(110) 



(111) 



It includes two types of local operators A v and Bf. 
The operator A v at vertex v defined by 



K = \G\~ X Yl L M> 

(l„£G 



(112) 



is an average of all local gauge transformations at v. This 
is the same as Eq. (| 108[) By checking that L v (h' v )A v = 
IGI" 1 E^ec Lv(KK) = IGI" 1 J2h v eG L M = A v , we 
see this is the projector that projects onto states that are 
invariant under the local gauge transformation L v (h' v ) at 
vertex v for any h' v e G. Therefore, A v prefers gauge 
symmetry at vertex v. While the gauge symmetry bro- 
ken states are allowed, it costs a energy of 1 to break 
the gauge symmetry. An important consequence of this 
gauge symmetry breaking is that a quantum number 
emerges at vertex v, and it is classified by the repre- 
sentations of the gauge group G. This quantum number 
identifies a quasiparticlc at vertex v. The group element 
g e represent the action on the states of the parallel trans- 
port of this emergent quasiparticlc along the edge e of the 
graph. 

The operator Bf on face / is defined via 



B t 




= 5, 



31 -32 -33 




(113) 



which is the same as Eq. ©. 

Here (715293 is the holonomy around the face /, and 
the delta function 5 a = 1 if the group element a equals 
the identity element e in G and otherwise. The delta 
function can be expanded in terms of characters, 



\G\~ 



E 

pGlrrep(G) 



dim p x p (g), 



where Irrep(G) is the set of all irreducible representations 
of G, dinip the dimension of the representation p, and x p 
the character of the p. Thus Bf is a projector that mea- 
sures whether the holonomy around the face / is trivial 
or not. 

Returning to the cases where a is in general nontrivial, 
our model may be viewed as the twisted version of 
Kitaev's QD model. In this interpretation, A 9 V is the ac- 
tion of the twisted gauge transformation at v, and the A v 
is the average of all local twisted gauge transformations. 
To make this interpretation precise, we need to study the 
algebra of all local operators, which is the main task of 
the next section. 
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IX. RELATION TO DIJKGRAAF— WITTEN 
TOPOLOGICAL GAUGE THEORY AND 
CONFORMAL FIELD THEORY 

In this section, wc dwell on the relation between our 
model, a lattice realization due to Dijkgraaf and Wittcn 
of topological Chern-Simons gauge theories, and confor- 
mal field theories. 

We begin with a quick review of the gist of the part 
of Dijkgraaf- Witten gauge theories that is relevant to 
our model. In Re^, Dijkgraaf and Witten established 
a correspondence between the three dimensional Chern- 
Simons gauge theories with a compact gauge group G and 
the two dimensional sigma models with Wess-Zumino 
interactions of the group G, in the sense that there is 
a natural map from the cohomology group H i (BG,7 t ), 
which classifies the Chern-Simons theories, and the 
group H 3 (G,Z), which classifies the Wess-Zumino inter- 
actions. The classifying space of the group G is denoted 
by BG. In general, the prescription of the topological ac- 
tion of a three dimensional Chern-Simons gauge theory 
is rather abstract; however, in view of that H 4 (BG, Z) is 
isomorphic to H i (BG 1 U(l)) when G is finite, Dijkgraaf 
and Witten constructed a concrete lattice realization of 
the topological action in the case of finite gauge groups. 
From now on in this section, we restrict the discussion to 
finite groups only. 

So, more precisely, consider a topological gauge the- 
ory defined in a three dimensional manifold M, with a 
finite gauge group G, the lattice realization is defined on 
a three-skeleton, i.e., a triangulation T, of M, with a 
group element of G living on each 1-simplex, which is 
oriented, of the triangulation (see Fig. The topo- 
logical partition function of such a lattice gauge theory 
reads 

Z(T(M),G) = -j-l[W(T i p, (114) 

' ' i 

where the product runs over all of the tetrahedra Tj 
in the triangulation T(M), and Ei is a sign, +1 or 
— 1, depending on whether the four vertices of the cor- 
responding tetrahedron are in a right-handed arrange- 
ment or left-handed arrangement. It is shown^l that 
the W(Tj) associated with each tetrahedron T, is a 3- 
cocycle over G. For example, for the tetrahedron in Fig. 
UW(T) = a(g,h,k). 

Note that a is an equivalence class, any two represen- 
tatives of the class are related by a 3-coboundary. If the 
manifold M is closed, the value of Z(T(M), G) does not 
depend on the choice of the representative of an equiv- 
alent class of 3-cocyclcs. The partition function is also 
invariant under the Pachner moves that connect two sim- 
plicial triangulations of M. Another remark is that the 
labeling of the vertices of the tetrahedron in Fig. [5] is 
fixed once for all. At this point, the partition function 
does not have tetrahedral symmetry. 

On closed manifolds, the partition function in Eq. 
(|114[) also has a gauge invariance. Consider the single 
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FIG. 6: A tetrahedron whose edges are graced with group 
elements; group multiplication rule applies to each of the four 
triangles. The corresponding 3-cocycle is a(g,h,k). 

tetrahedron in Fig. [S]as an example, the gauge transfor- 
mation that acts on vertex 1 transforms the topological 
action as follows. 

W(T) -> W'(T) = ^9,h)a(c 9h,k) w{n 
a(c,g, hk) 

where c £ G is the gauge parameter. This gauge transfor- 
mation can be understood topologically as in Fig. [7J The 



i 




FIG. 7: Gauge transformation acting on the vertex 1 in Fig. 
[6] a new vertex 1' is created at the barycenter, accompanied 
by a new group element c on the edge [l'l]. Also understood 
as the 1 — > 4 Pachner move. 

tetrahedron [1234] is associated with the original topo- 
logical action W(T). The gauge transformation acting 
on vertex 1 creates a new vertex 1' inside the tetrahe- 
dron (can be thought as being at the barycenter) and 
thus created four new tetrahedra, of which the tetrahe- 
dron [0'123] is associated with the new topological action 
W'(T) = a(cg, h, k). There are five tetrahedra all told in 
Fig. [7J associated with which the five 3-cocycles satisfy 
the 3-cocycle condition and thus lead to Eq. ()115p . On 
closed manifolds, the topological action is invariant under 
the gauge transformation Eq. pi5[) because the factors 
on the RHS of Eq. (19) can be canceled by those pro- 
duced by the gauge transformation on the neighbouring 
tetrahedra. Topologically, the gauge transformation be- 
haves like a 1 — > 4 Pachner move that splits a tetrahedron 
at it barycenter into four tetrahedra. Such a Pachner 
move can be visualized only in four-dimension, whereas 
Fig. [TJis the three-dimensional projection of a 4-simplcx 
whose five boundary 3-simplices are the five tetrahedra 
in the figure. 

To gain a deeper understanding of the gauge transfor- 
mation Eq. (|115|) . let us rewrite the equation in terms of 
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3-cocycles only as follows. 
ct(g, h, k) — ¥ a (g, h, k) 



a(cg, h, k) = 7(0, g, h, k)a(g, h, k), 



(116) 



where we define 



j(c,g,h,k) 



a(c,g, h)a(c,gh, k) 
a(c,g, hk) 



In general, however, the new object a'(g, h, k) is not a 3- 
cocycle any longer because one can check that it does not 
meet the 3-cocycle condition Eq. J3J). Nonetheless, that 
a! is not a 3-cocyclc makes it possible to choose a conve- 
nient gauge such that the prescription of the topological 
partition function becomes simpler. Indeed, according to 
Re^ 7 -, depending on the divisibility of |G|, the following 
gauge of the 3-cocycles may be imposed. 



ol {g,g 1 ,h) = a'(g,h,h *) = 1. 



(117) 



Under this gauge, the ordering of the vertices of a tetrahe- 
dron is irrelevant; in other words, the topological action 
W'(T) acquires tetrahedral symmetry, in the sense that 
it is invariant under the change of the labeling of the ver- 
tices. Therefore, the 3-cocycle condition and the gauge 
in Eq. ||116|) are incompatible unless the 3-cocycle un- 
der consideration is equivalent to the trivial one, namely 
at [a ]. 

If the manifold M has a boundary (open or closed), 
however, the gauge transformation in Eq. Q115P ceases 
to apply because a boundary condition must be imposed 
on M, which fixes the boundary value of the embedding 
of M into the classifying space BG and hence forbids 
the Pachner moves that involve the boundary simplices. 
Effectively, there are now degrees of freedom that can- 
not be gauged away on the two-dimensional boundary 
dM of M. As such, the three-dimensional partition 
function turns out to be the wave function of the cor- 
responding boundary state at certain time. Since the 
Dijkgraaf-Witten Chern-Simons theory is a topological 
gauge theory, there is no nonvanishing Hamiltonian due 
to Legendre transform that can enable the notion of 
ground and excited states for the boundary states. In- 
stead, here we have only gauge-invariant and noninvari- 
ant states. In particular, however, if the manifold M 
has a closed boundary, e.g., a solid torus with its bound- 
ary a 2-torus, the remaining boundary states are auto- 
matically only those gauge-invariant ones. This can be 
understood by the standard technique of "gluing" and 
"sewing" if a topological quantum field theory, which in 
the current case gives the size of the Hilbert space on 
the boundary. As such, the Dijkgraaf-Witten partition 
function becomes the dimension of the Hilbert space of 
the gauge-invariant states on the closed boundary of M. 

At this point, one may ask if it is possible to construct a 
Hamiltonian on the closed boundary of M whose ground 
states happen to be the gauge-invariant boundary states 
of the Dijkgraaf-Witten theory in M. Yes, the Hamilto- 
nian of our TQD model turns out to be a positive answer 
to this question, as explained as follows. 



Staring at Fig. [T5] again as if it is a triangulation of a 
3— ball, then the the right-to-left projection of the four 
triangles comprising the boundary 2-sphere to the paper 
plane is the very graph in Fig. |2(b)| a basis graph of our 
twisted quantum double model. 40 When the correspond- 
ing state of the graph in Fig. |2(b)| is a ground state, the 
Bf operator is unity acting on any of the triangles, com- 
plying with the group multiplication rule on each of the 
triangles in Fig[SJ 

By comparing A% with c G G at vertex 3 on the graph 
in Fig. 2(b) to Eq. 1115) one may see the vertex operator 



A 9 V in the Hamiltonian of our model is formally identi- 
cal to the gauge transformation on the partition function 
in Dijkgraaf-Witten theory. In our model, the opera- 
tors A v evolve the states; however, we have A v = 1 on 
the ground states, which implies that they are invariant 
under the Dijkgraaf-Witten gauge transformation in Eq. 
(| 1 1 5 jl . If we embed the graph of the TQD model on a 
torus or any surface homeomorphic to it, say, d X x S , 
where X is homeomorphic to a disk, we soon see that the 
ground states of the TQD model are the gauge-invariant 
boundary states of the Dijkraaf-Witten theory defined 
on M = X X S 1 . Therefore, we can conclude that 



->cs 



(T(X x S 1 ),^) = GSD TQD (T(dX x S 



*))■ 



The results obtained in Sections IIVI and [V] fall into this 
latter case, which can be verified by comparing them to 
the corresponding results in Re f 27 ' 35 . 

Now that we have gone through the logic of this corre- 
spondence, we can claim that our twisted quantum dou- 
ble models may indeed be regarded as a valid Hamilto- 
nian extension of the Dijkgraaf-Witten discrete Chern- 
Simons theories. Although the discussion so far is re- 
stricted to 3— dimension, the correspondence described 
above may be readily generalizable to higher dimensions. 

In light of this correspondence, the gauge transforma- 
tion in Eq. Q116p implies that all twisted quantum double 
models but the one defined by [a ] do not have tetrahe- 
dral symmetry. 

On the other hand, the CS theories in a 3-manifold 
M also correspond to other two-dimensional theories on 
the boundary dM, namely the rational conformal field 
theories (RCFT), such as Wess-Zumino-Witten (WZW) 
models. The salient point of this correspondence is that 
the CS Hilbert space on dM is isomorphic by a canoni- 
cal identification to the space of holomorphic conformal 
blocks of the RCFT on dM, while the CS wave-function 
reproduces the fusion algebra of the holomorphic sector 
of the RCFT. Note that this correspondence is level by 
level, in the sense that the CS theory and the WZW 
interaction of the corresponding RCFT are both at the 
same level, say, k. 



'Note that in this case, there are four triangles in Fig. |2(b)| 
including the triangle [124] in the back. 
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What follows naturally is a correspondence between 
our twisted quantum double models with a finite group 
G and a type of RCFTs, namely the CFTs as orbifolds by 
the group G of a holomorphic CFT. The modular data of 
a G-orbifold is twisted by a 3-cocycle over G. This cor- 
respondence is also level by level, in the following sense. 
The third cohomology group over G, H 3 (G,U(1)), is a 
discrete group; in particular it is Z if G is a compact Lie 
group. Hence, one can label the equivalence classes as the 
elements of H 3 (G,U(1)) by integers, say, [k]. The situa- 
tion is similar when G is finite, in which case H 3 (G, U(l)) 
also becomes finite, indicating that there are only finite 
number of levels available. Hence, a twisted quantum 
double model defined by a 3-cocycle a G [k] corresponds 
to a G-orbifold twisted by the same a, which has a 
twisted modular data at level k too. As such, when 
a — a° G [0] , the twisted quantum double model is actu- 
ally untwisted, which is equivalent to Kitaev model, and 
thus corresponds to the usual untwisted G-orbifold. 

Let Hc, a be a twisted quantum double model on a 
torus and Cg.u a twisted toroidal orbifold of a holomor- 
phic CFT C. One can check that they correspond to each 
other in the respects tabulated row by row as follows. 





Cc,a 


States: \g,h) 


conformal blocks: hd 

g 


Ground State: | A, fj.) 


l-loop characters: 


Ground state degeneracy 


Number of primary fields 


S & T matrices in \ A, y) basis 


S & T matrices in basis 



TABLE IV: Correspondence between a TQD and a twisted 
G-orbifold. 



A few remarks on the table are in order. The equality 
between the GSD of the TQD model and the number of 
primary fields of the corresponding orbifold is not surpris- 
ing, as each primary field is associated with a highest- 
weight vector of an irreducible representation of the Vira- 
soro algebra that is annihilated by the positive modes of 
the algebra and thus can be thought as a "ground state". 
Since each l-loop character counts all descendants of a 
primary field, including the primary field itself, it natu- 
rally corresponds to a unique ground state of the TQD 
model. Like the TQD ground states \A, fj,}, The l-loop 
characters form a orthonormal basis, in which the T 
operator is diagonal. The conformal blocks hO projects 

onto this basis as 

which is precisely how a TQD state \g, h) projects onto 
the ground state basis, as in Eq. (jBBjl . 

Although we have been talking about the fractional 
topological numbers and statistics of the quasiparticles 
of our model, we do not have in hand the operators 
that can create or annihilate these quasiparticles, nor do 



we know the exact wave functions of these quasiparti- 
cles. Nevertheless, as a ramification of the correspon- 
dence with the orbifold CFTs, that we can study the 
topological numbers and statistics of the possible quasi- 
particles of our model by using only the modular matrices 
can be expected. This ramification is further propped by 
a similar correspondence between the (2+1) -dimensional 
Hamiltonian formulation of fractional quantum Hall ef- 
fect (FQHE) systems and two-dimensional RCFT, which 
maps the holomorphic wave functions of the quasiparti- 
cles of the FQHE system to the conformal blocks of the 
CFT. 



X. RELATIONS TO LEVIN-WEN MODELS 

In this section, we discuss the relation between Levin- 
Wen models and our TQD models. In particular, we 
demonstrate a duality map of a class of Levin- Wen mod- 
els into certain TQD models. 

To begin with, let us briefly review Levin-Wen Mod- 
els. Levin- Wen models, also known as string-net models, 
were proposed to generate the ground states that exhibit 
the phenomenon of string-net condensation as a physi- 
cal mechanism for the time reversal invariant topological 
phases. They are believed to be a Hamiltonian formu- 
lation of the Turaev-Viro topological field theories^—, 
analogous to that our TQD models arc a Hamiltonian 
extension of topological Chern-Simons theories, as be- 
labored in the previous section. Levin- Wen models are 
usually defined on the honeycomb lattice. 

String degrees of freedom reside on the edges of the 
honeycomb lattice, each link of which is graced with one 
of N + 1 string types. In the most general setting, the 
N + 1 string types form a finite set / equipped with a 
duality map * : I —> I such that j** = j for all j G 
J. These abstract string types are usually considered to 
label the irreducible representations of certain group or 
algebra (e.g., a quantum group). 

A Levin- Wen model is specified by a triple of the input 
data {d, 8, F}. Quantum numbers dj are called quantum 
dimensions and are complex numbers associated with the 
group elements j G G, satisfying dj = dj* . There are in 
principle two ways of setting up the fusion rules 5. First, 
one can let the tensor product rules of the irreducible 
represcnations labeled by the string types as the fusion 
rules. But we do not consider this case here. Second, 
which is the case to be discussed in this section, one can 
use the product rule of certain group G as the fusion 
rule, which is in fact a Kronecker (^-function associated 
with each triple of string types {i,j, k} respectively on 
the three links meeting at a vertex, such that 8^ equals 
1 if the group multiplication ijk is the identity clement 
e, and otherwise. The quantum dimensions and fusion 
rules must satisfy 

didj = Sijk*dk- (H8) 
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The dual string type j* = j _1 can also represent the cor- 
responding inverse group clement in G. Fig. [5] illustrates 
the fusion rule on the honeycomb lattice, by showing just 
one vertex. A link a of the lattice is graced with a string 
type s a and endowed with an orientation, specified by an 
arrow. Such a string of type s a can also be represented 
by a flipped arrow, but with the conjugate string type 

S* a = (Sa)- 1 . 





FIG. 8: A string-net vertex. 

The Qj symbol F arc complex numbers that obey the 
following self-consistency conditions 



thus the +1 eigenstates of A v and B p , which are known 
as the string-net condensed states. 

We emphasize that in Levin-Wen models, the plaquc- 
tte operator B p is identically zero outside of the subspace 

of A v = 1 because the F-symbols are automatically zero. 
This is in contrast to our TQD models in which the vertex 
operators are well-defined and nontrivial outside T-L Bf=1 . 

The honeycomb lattice has as its dual lattice the trian- 
gular lattice, which is a regular case of the triangle graph 
r on which our TQD models are defined. We further no- 
tice that the way we enumerate the vertices and assign 
group elements on the edges of T does induce orienta- 
tions on the links of the dual honeycomb lattice, as seen 
in Eq. (|122[) . in which only the dual honeycomb plaquc- 
ttes are shown. Bear in mind that the enumerations of 
the vertices of T now label the plaqucttcs. 



r kln 



pmij v m v n 
nk*l* VjVi 



pklm* 
r ijn" 



r lkn* 



r l*k*n 



rpmlq pjip pjs'n 
kp* n r mns* Ikr* 



pinlq pl* m *i* 
kp*n pk*n 



F 3W pr 

q* kr* mis* ' 



(119) 

where Vj = ^fdj (and v e = 1 for the identity group ele- 
ment e). The first line is a symmetry over the indices of 
F, where the last equality is meant for the Hamiltonian 
to be Hermitian, the second line the pentagon identity, 
and the last line the orthogonality condition. We re- 
mark here that it is the solution of the _F-symbols and 
the quantum dimensions to Eq. (|119[) that dictates for 
which the abstract string types label the irreducible rep- 
resentations. Often, the string types turn out to label 
the irreducible representations of a rather complicated 
algebra 21 although the group G that supplies the fusion 
rules is very simple. See Section [VII Bl for an example. 
The usual string net Hamiltonian takes the form 



h = — a v — b p , 



(120) 



where A, 



= S. 



* j it- 



is the vertex operator de- 



fined at each string-net vertex, and B p = ^2^ = Qa s B p} 



with a s = d s /D and D = J^iLo^f: is the "magnetic- 
flux" operator defined for each hexagonal plaquette of 
the string-net lattice. Each operator B p in B p acts on a 
plaquette as follows^. 



Bi 




a-ie'; 




(121) 



where e' a = e a s. 

The vertex operators A v are projectors. The parame- 
ter a s ensures that the operators B p are also projectors. 
It can be shown that {A v , B p \ \/v, p} is a set of commuting 
operators, whose common eigenstates span the Hilbcrt 
space of the model. The ground states of the model are 



''4 



('1 



V-2 



[VlV 2 




(122) 



The pair [vii^] indicates that the plaquette v\ is on the 
left of their common edge, while V2 is on the right. Hence, 
the group element on the edge can be denoted as gi VlV2 ] 
and obviously satisfies g[ Vl v 2 ] = 9[vl Vl ]- This rule can be 
applied to the entire graph. 

These observations imply that there may exist a kind of 
duality between the concerned type of Levin- Wen mod- 
els and certain TQD models, in fact an inclusion of the 
former into the latter, as we now explore. 

We claim that any _F-symbol that solves Eq. (|119p can 
be mapped to a 3-cocycle a (up to a 3-coboundary) that 
defines a TQD model. That said, a Levin-Wen model 
with such an F can be identified with a TQD model 
with the corresponding a. Indeed, a tensor F has three 
independent indices and can be expressed as 



r kln 



a(i,j, k)6 m ^ifi-iS n! jk5mjk)-i 



(123) 



where a 6 C x = C \ {0} is a function of i, j, k g G. The 
self-consistency conditions in Eq. (|119|) become 



VijVjk 
VjV ijk 



a(i,j,k) =a((ij) 1 ,i,jk) 
= a (k, (ijk)~ l ,i) 

a(l,k,j) a(m,lk,j) 



a{ml, k,j) a(m, I, kj) 
^a(m, I, fc)a(/ _1 , m , mlk) 



a(m, I, k) = 1, 
1. 



(124) 

respectively. The pentagon identity in the second line is 
readily the 3-cocycle condition of this a that turns out to 
be a 3-cocycle in H 3 (G, C x ). It does not harm to assume 
that a is [/(l)-valued in H 3 (G,U{1)) ^ iT 3 (G,C x ). 
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Then we see that the orthogonality condition is identified 
with one equality in the symmetry condition. 

The quantum dimension dj then takes the definition 

d i = a(r 1 ,i,i- 1 ), (125) 

as is verified by setting j — i^ 1 and k — i in the symmetry 
condition and by the 3-cocycle condition. 

The conditions didk = d^ and di = dj-i = ±1 are 
immediate consequences of the symmetry condition in 
Eq. (|124p . By setting j to the identity element e in the 
first equality of the symmetry condition (|124l) and by 
applying the 3-cocycle condition 5a ((ij) -1 , i, j, k) = 1, 
we obtain 

ViVk a(e,e,k) 
Vih a[i 1 ,i,k)a(i L ,i,e) 

Since a(z _1 ,z,fc) 2 = 1 from the symmetry condition for 
all i,k G G, taking a square of the above equation yields 
didk = djfc. The condition di = d^-i = ±1 is due to 
a(i , i, i) 2 = 1 and djdj-i = d e = 1. 

Substituting the a in Eq. (|123|l into the in Eq. 

(|12ip . we find that A*|b /= i can be identified with d s i?p 
for all p up to a common unitary transformation that de- 
pends on the enumeration of the vertices and the choice 
of {vi}, where A"\B f =i is the TQD model vertex operator 
acting on the vertex p dual to the plaquette p on the hon- 
eycomb lattice, with the action restricted to the subspacc 
H Bf=1 . Notice that the enumeration dependence of the 
transformation is related to the enumeration dependence 
of the definition of A s p . 

Therefore, we infer that the Levin- Wen model with the 
fusion rule 5ijk determined by the group multiplication 
law of G and the F-symbols meeting Eq. (|119[) can be 
identified as the TQD model restricted to H Bf=1 that is 
defined by the (not necessarily normalized) [/(l)-valucd 
3-cocycle a that satisfies 

a(i,j,k) = a jk) 13 ° k = a (k, (ijk)~ l ,i) 

= a(j,i,ijk- 1 )=a(j- 1 ,i- 1 ,ijk)~ 1 (127) 

where Vi — s/dl and di = a(i,i _1 ,i). 

Nevertheless, not all 3-cocycles satisfy the conditions 
(|127(1 . in which case the TQD models do not correspond 
to any Levin- Wen models. For example, with Z2 fusion 
rule, there are two Levin-Wen models, namely the toric 
code model and the double semion model, which can be 
identified with the TQD models respectively with the [a ] 
and [a 1 ] given in Eq. (|89[) . which correspond to di = 1 
and di = — 1 respectively, according to Eq. Q125P . There 
is, however, only one Levin- Wen model with Z3 fusion 
rule that has a dual TQD model — the one with the [a ] 
given in Eq. (f9U|) . Note that here, we do not distin- 
guish the TQD models defined by equivalent 3-cocycles 
because they describe the same topological phases, as ex- 
plained in Section III CI Details of the TQD models with 
Z2 and Z3 are found in Section IVlI Bl 



The study in this section partially answers the question 
when and how Levin- Wen models can be characterized 
by group cohomology, which was raised in Ref^ 2 -. 



XI. DISCUSSIONS AND OUTLOOK 

In this very last section, we shall summarize our major 
results along with discussions on a few questions bonded 
to these results that are yet not fully answered in this 
paper but deserve future exploration. 

First of all, we fabricated a new model — the Twisted 
Quantum Double model — of 2d topological phases by a 
3-cocycle [a] e H 3 (G,U(1)) of a finite group G on a 
graph composed of triangles, each edge of which is dec- 
orated by an clement of G. This model constitutes a 
very rich class of topological phases, which are otherwise 
missing in some other models, such as the Kitaev model 
and Levin- Wen model. The topological properties of the 
TQD model are reflected in the topological numbers — 
GSD, topological spin, etc — associated with the topolog- 
ical observables of the model. 

We further classified these topological numbers, which 
cither directly depend on the defining 3-cocycle of the 
model or indirectly via a twisted 2-cocycle determined by 
the 3-cocycle. Two TQD models defined by two equiv- 
alent 3-cocycles are shown to bear the same topological 
phase. We thus expect that the classification of the topo- 
logical numbers does the job as well for the topological 
phases described by the TQD models. We expect but do 
not affirm this yet because we have not been able to ex- 
plicitly prove that two inequivalent 3-cocycles never yield 
the same topological phase. This and detailed studies of 
the topological phases certainly calls for more efforts in 
future works. 

Second, our TQD model appears to be certain gener- 
alization of the Kitaev model in the following sense. A 
TQD model is precisely a Kitaev model when its defin- 
ing 3-cocycle is trivial. In this situation, the Hamilto- 
nian consists of local gauge transformations and local 
flux projections. As a collective effect, the ground states 
are classified by the irreducible representations of quan- 
tum double of the finite group G, and this is expected to 
be true also for the quasiparticlc excitations. When the 
defining 3-cocyclc is non-trivial, the Hamiltonian can be 
viewed as consisting of local twisted gauge transforma- 
tions and local flux projections. Similarly the ground 
states are classified by the irreducible representations of 
the twisted quantum double of the finite group G. 

Third, we relate our TQD models to the Dijkgraaf- 
Wittcn (DW) topological Chcrn-Simons theories, by 
viewing ours as a Hamiltonian extension of the latter. In 
fact, we have shown that the GSD of a TQD model de- 
fined by some 3-cocycle on the boundary of a 3-manifold 
coincides with the partition function of the DW topolog- 
ical Chcrn-Simons theory in the bulk, whose topological 
action is given by the same 3-cocycle. 

This connection motivates a correspondence between 
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our TQD models on a torus and the RCFTs that are 
the toric orbifolds by a finite group of a holomorphic 
CFT and are twisted by nontrivial 3-cocycles. This 
correspondence identifies the ground states, the GSD, 
and the modular matrices of a TQD model, respectively, 
with the holomorphic characters, the number of primary 
fields, and also the modular matrices of the correspond- 
ing RCFT. Provided with the description of fractional 
quantum Hall effect by CFT, we are encouraged to ex- 
pect that the statistical and topological properties of the 
quasiparticle excitations and hence the topological phase 
of a TQD model can be investigated in terms of the mod- 
ular matrices of the model. 

Third, to echo the fact that our model is partly moti- 
vated by the Levin-Wen model, we studied the relation 
between TQD models and the type of Levin- Wen models 
where the fusion rules coincide with the multiplication 
laws of finite groups: we demonstrated that each such 
Levin- Wen model on a graph can be directly translated 
to a TQD model of the type on the dual graph. The 
reverse is not true, however, indicating that TQD mod- 
els embodies more topological phases than this type of 
Levin- Wen models. In our study of this in Section [Xl we 
adopted simply the original settings of the Levin-Wen 
modelii. 

Furthermore, as pointed out in SectionlXland in Re^, 
the fusion rules in a Levin- Wen model can in principle 
be identified with the tensor product rules of the irre- 
ducible representations of certain group or algebra. If 
this is the case, the pentagon identity in Eq. (|119| of 
the _F-symbols contains a summation over the index n 
on the LHS due to the summation that would appear 
in the fusion rules, which comprises the interpretation of 
the F-symbols as 3-cocycles and the pentagon identity 
as the corresponding 3-cocycle condition. Then clearly, 
this type of Levin-Wen models, apart from the special 
cases where the representations are restricted so as to re- 
move the summation, cannot be dual to our TQD mod- 
els. This type of Levin-Wen models are believed to be 
classified in terms of tensor categories. This distinction 
between the two types of Levin- Wen models is related 
to the question when and how Levin-Wen model can 
be classified by group cohomology, which is raised in 
Ref22 that is inspired by the duality between certain SPT 
phases and long range entangled topological phases de- 
scribed by certain Levin- Wen model a 22 ' 427 . The duality 
found between our TQD models and Levin- Wen models 
then partially answered this question. 

Above all, a main purpose of this paper is to reveal 
the topological properties of the ground states. We pro- 
pose the following topological properties of the elemen- 
tary excitations. The number of quasiparticle species in 
the elementary excitations is equal to the GSD on a torus. 
Moreover, the topological charge that identifies the quasi- 
particles are classified by the twisted quantum double of 
the finite group G, and The S and T statistical matrices 
are the same as the modular S and T matrices derived 
from the topological observable in the ground states on 



a torus. Work is in progress in this direction. 
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Appendix A: Algebra of local operators 

In this appendix, we show that the local operators in 
the Hamiltonian Eq. ([5]) forms the following algebra. 

(i) [B f ,,B f ]=0, [B f ,A«]=0, 

(ii) [Al,A h w ] = Qiiv^w, 

(iii) Al,A h v =A% h , (Al) 

where in (iii), [v'v] = g is understood. The equality (iii) 
in the above implies that A v = J2 g -^2/1^1 i s a projector. 

(i). That [Bf,Bf] = follows immediately from the def- 
inition of Bf in Eq. Since A v affects only the group 
elements on the boundary of the vertex v, [Bf,A^] = 
holds obviously when the vertex v is not on the bound- 
ary of /. In the case where v is right on the boundary of 
v, without loss of generality, let us consider two actions, 
A 9 V2 Bf and BfA^ 2 , on the basis vector 




We have 




(A2) 







[v' 2 v 3 ] 






' a([«i^], 


[v' 2 v 2 


, [V2V3]) 




3[viv 2 ] 


[v 2 v 3 \ 


[V3V1] 






[v' 2 v 2 


, [V2V3}) 






[v 2 v 3 ] 


[V3V1] 



a{[viv' 2 ), [v' 2 v 2 ], [v 2 v 3 ]) 
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(A3) 



where [v' 2 V2] — g is understood, the dots . . . collects 
all other factors irrelevant and thus omitted, and the 
second equality follows from applying the chain rule 
[vrv' 2 ] ■ [v' 2 v 3 ] = [v l v' 2 ]-[v' 2 v 2 }-\v2v' 2 ]-[v' 2 v i ] = [wit^Hi^]- 
Hence, we conclude that [Bf,A%] = holds for all 

f,ver. 

(ii) . It is clear by the definition of B v that if v\ and v 2 
are not connected by any edge, [A^Aj-J = is true. 
We then need only to check the case where v\ and v 2 are 
neighboring to each other. Let us first check the following 
action of A 9 Vi A^ 2 on a relevant basis vector. 



A9 A h 



--{a([viv 2 }, [v' 2 v 2 ], [v 2 V4,]) 1 



:([V1V 2 }, [v' 2 V2], [V2V 3 ]) -)A 9 V 



-(a ([viv 2 ], [v' 2 v 2 ], [v 2 Vi]) 1 a {[viV 2 ], [v 2 v 2 ], [v 2 v 3 ]) 
(a([v[vi], [viV 2 ], [v' 2 v 4 ]) 



:(Kui], [viv 2 ], [v' 2 v 3 }y 



(A4) 



with [v' 2 v 2 ] = h, [v{vi] = g. Only those 3-cocycles corre- 
sponding to the two common boundary vertices are writ- 
ten down. 

By using twice the 3-cocycle condition in Q , we have 

"(K^l: W 2 v 2 \, [v2Vi]) a ([v'^} • [uiUj], [v' 2 v 2 ], [v^^y 1 
a ([v[vi], [vtv^] ■ [v' 2 v 2 ], [v 2 v 4 ]) 

a ([v[vi], [viv' 2 ], [v' 2 v 2 ] ■ [v 2 v 4 ]) 1 x 
a([v[vi], [viv' 2 ], [v' 2 v 2 ]) = 1 



(A5) 



and 



«(K u 2L ["2^2], [V2V3]) 

a (K«i] • {viv' 2 }, [v' 2 v 2 ], [v^])' 1 
a ([v[vi], [vW^ ■ [v' 2 v 2 ], [v 2 v 3 ]) 

a ([v[vx], [viv 2 ], [v' 2 v 2 ] ■ [v 2 v 3 })~ 
a([v[vi], [vW 2 \, [v' 2 v 2 ]) = 1 



(A6) 



of A 9 V1 A% 2 is the same as A^ 2 A 9 ]1 : 



A h A 9 



a([w^i], [viv 2 ], [V2V4}) 



a([v[vi], [viv 2 ] 7 [u2f4])a(Kwi], [viv 2 ], [^2^3])" 
"(K^], W 2 v 2 ], [v 2 va]) 1 



<X{[V' 1 V 2 ],[V 2 V 2 ],[V2V3]) ■■ 




(A7) 



Notice that the chain rule in Eq. (|10p guarantees that 
each group element indexed by the same pair of enumer- 
ations is the same in the above evaluations. Therefore 
we arrive at A^A^ = A^A^. 

(iii). We show A 9 v ,A h v = A g v h , where v is assumed to 
become v' after the action of A% with [v'v] = h, while 
v is turned to be v" with [v"v] = g ■ h. Wc begin with 
the action of A^A^ on the vertex z; 2 of the basis vector 
as follows. 



A a ,A h 

v' V2 




=a([viv' 2 ], [v' 2 v 2 ], [v 2 v 3 ]) 

a ([ v lV 2 }, [V 2 V 2], {V2V4W 1 

a ([ V 2 V 2], [V2V 3 ], [vsViiy 1 A 9 V , 



{a{[viv 2 }, [v' 2 v 2 ], [v 2 v 3 ])a{[viv 2 ], [v' 2 v 2 ], [v 2 Vi]) 1 
a([v2V 2 ], [v' 2 v 3 ], [v-iV^y 1 ^ 
a [v 2 W 2 ], [v' 2 v 3 })a([ Vl v' 2 l [v' 2 %], [v'My 1 




"(K^L [«2 U 3], [U3V4]) 1 




(A8) 



together with the chain rule Eq. (|10p. we find the action 



Using three times the 3-cocycle condition J3| , 

u{[viV 2 ), [v 2 v 2 ], [v2V 3 })a([v 1 v 2 ], [v 2 v 2 ], [v' 2 v 3 ]) 
=a([v 1 v 2 ), [v 2 v 2 ], [v' 2 v 2 ])a{[viv' 2 % [v' 2 v 2 ], [V2V3]) 
a([v' 2 W 2 },[v 2 v 2 ],[v 2 v 3 }) (A9) 
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a ([ v lV 2 ], K^], [ W 2 W 4]) a ([v"v' 2 ], [v' 2 V 2 }, [V2V4}) 

a(K^],K« 2 ],[« 2V4 ]) _1 (A10) 



(|B2|) we can neglect the part of T out of the scope of the 
action of Tj and hence the A v acting on this part. Wc 
now pursue the proof in the following equations respec- 
tively for T\, T 2 , and T3. In the sequel, [ijk] denotes a 
triangle whose vertices are i, j, and k counter-clockwise, 
and "• • • " represents all the irrelevant factors, which are 
thus omitted. 



a{[v' 2 v 2 ],[v2V3],[v 3 V4}) 1 a([v 2 v 2 },[v 2 v 3 ],[v 3 v 4 ]) 1 

= a (K^l: [ v 2 v 2\i [V2V3})' 1 a {[V 2 V 2 ], [v' 2 V 2 ], [v 2 V 4 ]) 

a (K^], [V2V3}, [V 3 V 4 ]) (All) 

we obtain 

a ([vW^j, [v 2 v 2 ], {v 2 v 3 })a([viv 2 }, [v 2 v 2 ], [i^])" 1 



a ([v 2 v 2 ], [v2V 3 ], [v 3 v 4 ]) 




(A12) 



According to the chain rule in Eq. (fTU|) . we have 
v' 2 v 2 ]. The above action is identified 



[v 2 v 2 \ = [v'^v' 2 \ ■ 

as action of A^' o h . Therefore we conclude that 



49 A h — A9 h 
A v' A v — A v 

Though the above proof is done on a triangle plaquettc, 
the general proof on a plaquette of any other shape is 
straightforward. 



Appendix B: Mutations are Unitary Symmetry 
transformations 

1. Symmetry 



i=i 



\G\ 



- y 

3 4 L^i 



■a([l'l],[12],[24]) 



[1'1],[2'2], 
[3'3],[4'4]eG 



x«([l'2'],[2'2],[24])- 1 a([2'2],[23],[34]) 
xa([2'3'],[3'3],[34])- 1 «([l'2'],[2'4'],[4'4]) 

4' 

x a([2'3'],[3'4'],[4'4])Ti 



1 



\G\ A 



E 



•a([l'l],[12],[24]) 



[1'3'],[1'1],[2'2], 
[3'3],[4'4]SG 



xa([l'2'],[2'2],[24])- 1 a([2'2],[23],[34]) 
xa([2'3'],[3'4'],[4'4])a([l'2'],[2'3'],[3'4']) 

4' 




(B3) 



When restricted to ground states in C \ we 

can impose the following chain rule on all triangles in T: 



[ViVj] = [ViVj] ■ [VjV k ], 



(Bl) 



where Uj, Vj, and Vk are the three vertices of any triangle. 
Since the mutation operators are defined on the subspace 
np =l for all T, to show that TjP r = P r ,T t , where T = 
Tj(T), we can neglect the face operators in the projectors 
Pr and Pp' defined in Eq. (jT^J)- We thus need to show 
that 



Ti n k = n a ^ 



(B2) 



holds for all mutation operators Tj, i — 1,2,3, and any 
state in Hp /_1 on any graph T. Since a Tj acts on at 
most four three triangles and does not affect any other 
triangle in the same graph, Tj certainly commutes with 
any A v at any vertex v that does not lie on the bound- 
ary of the triangles on which the Tj acts. Hence, in Eq. 



([124], [234]) 



i=l 



: a ([12], [23], [34]) J] 4 
[13] eG 




E 



•a([l'l],[12],[23]) 



[13],[1'1],[2'2], 
[3'3],[4'4]eG 



xa([l'l],[13],[34])a([l'2'],[2'2],[23])- 1 
xa([l'2'],[2'3'],[3'3])a([l'3'],[3'3],[34]) 
xa([l'3'],[3'4'],[4'4])a([12],[23],[34]) 

4' 



(B4) 
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It is straightforward to show that the RHS of Eq. rules [l'l] ■ [13] = [1'3] and [1'3'] ■ [3'3] = [1'3], and by 
()B3ll is equal to that of Eq. (|B4j) by knowing that applying the following four 3-cocycle conditions in order. 



E 



[1'1],[13],[3'3] 



E 



[1'1],[1'3'],[3'3] 



because of the chain 



a([l'2'], [2'4'], [4'4])a([2'3'], [3'4'], [4'4]) = a([l'3'], [3'4'], [4'4]) a ([1'2'], [2'3'], [3'4]) a ([1'2'], [2'3'], [3'4'D" 1 , 

a ([l'l], [12], [24]) = a ([12], [23], [34]) a ([1'2], [23], [34])" 1 a ([l'l], [13], [34]) a ([l'l], [12], [23]) , 

a ([2'2], [23], [34]) a ([1'2], [23], [34])" 1 = a ([l'2'] ; [2'3], [34])" 1 a ([1'2'], [2'2], [24]) a ([1'2'], [2'2], [23]) , 

a ([2'3'], [3'3], [34])" 1 a ([1'2'], [2'3'], [3'4]) a ([1'2'], [2'3], [34]) _1 = a ([1'3'], [3'3], [34]) _1 a ([1'2'], [2'3'], [3'3]) . 



Thus, Eq. dB2j) holds for Ti. 

The case of T 2 is a bit trickier. Let us write down how 



the LHS and RHS of Eq. (|B2[I act on a state as follows 



J 



T 2 ([2'3'4'])n^ 



i=i 




pa ° t 23 i' t 34 i) ° (t 2 ' 3 ']' m t^r 1 ° (t 2 ' 3 '^ t 3 ' 4 ']' t 4 ' 4 ]) ^ 



[2'2],[3'3], 
[4'4]eG 




2 a([2'2],[23],[34])a([2'3'],[3'3],[34])- 



[2'2],[3'3],[4'4], 
[12'],[13'],[14']eG 



x a ( [2'3'] , [3'4'] , [4' 4] ) a ( [12'] , [2'3'] , [3'4' ] 




(B5) 



l[A t T 2 ([234}) 




= ]T « ([12], [23], [34])-^ A 

[12], [13], i=l 
[14] GG 




|G| 



^ £ ([12], [23], [34]) a ([l'l], [13], [34]) a ([l'l], [12], [24])" 1 a ([l'l], [12], [23]) 



[12], [13], [14], [l'l], 
[2'2],[3'3],[4'4]eG 



x a ([1'2'], [2'2], [23])- 1 a ([1'2'], [2'2], [24]) a ([1'2'], [2'3'], [3'3]) 



x a ([1'3'], [3'3], [34])" 1 a ([1'2'], [2'4'], [^J)" 1 a ([l'3'], [3'4], [44]) 




E 



([2'2],[23],[34])a([2'3'],[3'3],[34])- 



[2'2],[3'3],[4'4], 
[1'2'],[1'3'],[1'4']GG 
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x a ([2'3'], [3'4'], [4'4]) a ([1'2'], [2'3'], [3'4']) 




(B6) 



where the last equality is obtained by first plugging into 
the second row the following four 3-cocycle conditions 

q([l'l],[13],[34])q([l'l],[12],[23]) _ q([l'2],[23],[34]) 
q([l'l],[12],[24]) - q([12],[23],[34]) 

q([l'2 , ],[2'2],[24]) _ q([2'2],[23] [34])a([l'2'] [2'3],[34]) 
q([l'2'j,[2'2j,[23j) q([l'2j,[23J,[34j) 

q([l'2'],[2'3'],[3'3]) _ q([l'2'],[2'3'],[3'4]) 
q([l'3'],[3'3],[34]) ~ q([l'2'] ,[2'3] ,[34])q([2'3'] ,[3'3],[34]) 

q([l'3'],[3'4'],[4'4]) _ q([2'3'],[3'4'],[4'4])q([l'2'],[2'3'],[3'4']) 
a([l'2'],[2'4'],[4'4]) ~ a([l'2'],[2'3'],[3'4]) 

and then by applying to the summations these chain 
rules: 

[l'l] • [12] • [2'2] _1 = [1'2'], 
[l'l] ■ [13] ■ [3'3]- x = [1'3'], 
[l'l] • [14] ■ [4'4]-! = [1'4'], 



which arc guaranteed by the restriction of Bf = 1 and the 
properties of the vertex operators Ai. The summations 
over [2'2], [3'3], and [4' 4] can be replaced by those over 
[1'2'], [1'3'], and [1'4'] respectively. Since [l'l] does not 
appear in the 3-cocycles any more, Sri'i] contributes a 
factor a |G|. Clearly, the vertex enumerated by 1' on the 
RHS of the last equality of Eq. (|B6[) is now dummy and 
thus can be re-enumerated by 1 without altering its order 
relative the enumerations of the other three vertices. As 
such, we can see that the RHS of Eq. jB6| and that of 
Eq. l|B5p are actually identical. That is, Eq. (|B2[) is true 
for T 2 . 



When it comes to the case of T3 , we have the action of 
the LHS and that of the RHS of Eq. (|B2) on the same 
state respectively being 



T 3 ([l'2'3'4'])nA, 




J2 ' ' ' a (I 1 ' 1 ]' I 12 !' I 24 ]) a ([l'l]. [12], [23])" 1 a ([1'2'], [2% [23]) 



[1'1],[2'2], 
[3'3],[4'4]eG 



x a ([1'2'], [2'2], [24])' 1 a ([2'2], [23], [34]) a ([1'2'], [2'3'], [3'3])" 



x a([2'3'], [3'3], [34])" 1 a([l'2'], [2'4'], [4'4]) a ([2'3'], [3'4'], [4'4])T 3 




\G\ 



L J2 ■ • • « ([l'l], [12], [24]) a ([l'l], [12], [23])- 1 a ([l'2'], [2'2], [23]) a Ql'2'], [2'2], [24])" 1 



[1'1],[2'2], 
[3'3],[4'4]GG 

x a ([2'2], [23], [34]) a ([l'2'], [2'3'], [S'S])" 1 a ([2'3'], [3'3], [34])" 1 
x a ([1'2'], [2'4'J, [4'4]) a ([2'3'], [3'4'], [4'4]) a ([1'2'J, [2'3'], [3'4']) 




(B7) 
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([12], [23], [34]) [] A 




(B8) 



= -^a ([12], [23], [34]) £ a ([l'l], [13], [34]) a ([l'3'], [3'3], [34])" 1 a ([1'3'], [3'4'], [4'4]) 



[1'1],[3'3], 
[4'4]eG 




r 



The RHS of Eq. (|B7)) and that of Eq. ([B8]) can be the corresponding 3-cocycles in Eq. |(B7|) . 
identified by applying the four 3-cocycle conditions to 



«([l'l],[12],[24])a([l'l],[12],[23])- 1 =a([12],[23],[34])a([l'l],[13],[34])«([l'2],[23],^^ 

a ([1'2'], [2'2], [23]) a ([l'2'], [2'2], [24])" 1 a ([2% [23], [34]) = a ([1'2], [23], [34]) a ([1'2'], [2'3], [34])" 1 

a ([1'2'], [2'3'], [S'3])- 1 a ([2'3'], [3'3], [34])" 1 = a ([1'3'], [3'3], [34])" 1 a ([1'2'], [2'3], [34]) a ([l'2'], [2'3'], [S'4])- 1 

a ([2'3'], [3'4'], [4'4]) a ([1'2'], [2'4'], [4'4]) = a ([l'3'], [3'4'], [4'4]) a ([l'2'], [2'3'], [3'4]) a ([1'2'], [2'3'], [3'4']) _1 . 

I 



(B9) 



These 3-cocycle conditions render the 3-cocycles resulted 
independent of the J2[2'2] m ^1- fEZj, which then con- 
tributes a factor of \G\. Therefore, we conclude that Eq. 
(|B2[) holds for all mutation operators and any state in 

B i — 1 

H r , which means that the mutation operators pre- 
serve the space of ground states. 



2. Unitarity 

Now we need to show that all mutation transforma- 
tions are unitary, i.e., they satisfy Eq. (j2~4")l . It is suffi- 
cient to show all generators T\, and T3 are unitary. 
We show them respectively in the following. 

We first demonstrate that T\ is unitary not only on 
the ground states but also so over the entire subspacc 
H Bf=1 . Let us consider the action of on a generic 
basis state as follows, in which only the relevant part of 
the graph is shown. 



[«i«s]GG 

yi E a (N" 2 ]' [^2^3], [V3V4]) 

[v 1 v 3 ]eG [v 2 v 4 yeG 



x a {[V1V2], [V2V3], [V3V4]) 1 



E E 




E 

[«iu 3 ]eG 




which deserves some explanation. In the second equality, 
the action of T\ on two triangles sharing a horizontal 
edge is understood from rotating by 7r/2 either way the 
action of T\ on two triangles with a vertical common 
edge as in the first equality because there is a global 
rotation symmetry on the graph. The inverse 3-cocycle 
a -1 follows the rule described below Eq. (f2"T]) . The group 
element [i^imJ' that is summed over in the second equality 
is brought by the action of the second T\, as the edge 
[^2^4] in the second equality is a new edge relative to 
the [W2W4] on the LHS of the equation, despite the same 
end vertices, and thus should be graced with a different 
group clement [1)2^4]', which is however constrained by 
Bf = 1 on the four triangles to be equal to [^2^4]. The 
last equality is also due to the restriction to the subspacc 
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H B s- x . As a result, T± =T X , such that we can define 
that t\ = T\ and infer that T\ is unitary over the entire 
H B f =1 . 

Next, we show that T3T2 = 1 on the entire subspace 
T~L Bf=1 . We consider the action of T 3 T 2 on a generic 
basis state as follows, in which only the relevant part of 
the graph is shown. 



T0T0 



=P°(1,2,3) a([12'],[2'2],[23]) 




= X] a (k v l}d v l v 2},[ v 2V3})T 3 

[viq],[v 2 q], 

[v 3 q\eG 

= Qf([g«i],[ui« 2 ],[«2«3]) 

[viq],[v 2 q],[v 3 q]eG 



x u([<l v i}, [viv 2 ], [V2V3]) 1 






where the inverse 3-cocycle a -1 follows from the rule 
introduced below Eq. ([231) . At this point, one may think 
that T3 is the inverse of T 2 on H Bf=1 . But this is not true 
because T 2 T 3 ^ 1 in general. Nevertheless, fortunately, 
as we now show, T2T3 = 1 on the ground states T~L°. 
Since T 2 T 3 P° = T 2 P°T 3 = P°T 2 T 3 on n°, we have 



r 2 r 3 P (l,2,3,4) 




=r 2 P (l,2,3,4)r : 



=a([12],[23],[34])T 2 P°(l,2,3) 



=a([12],[23],[34])P°(l,2,3)T 2 



=P°(1,2,3) ]T a ([12], [23], [34]) 

[12'],[2'3], 
[2'4]£G 



xa([12'],[2'3],[34])- 





[12'],[2'3], 
[2'4]eG 



xa 



([2'2],[23],[34])a([12'],[2'2],[24])- 1 




=P (1,2,3)|^ «([12'],[2'2],[23]) 

' ' [2'2]6G 

xa([2'2],[23] ) [34])a([12'] ) [2'2],[24])- 1 



=P°(1,2,3)A 




=P°(1,2,3,4) 




(BIO) 



where P°(l,2,3) and P°(l,2,3,4) are projectors acting 
on the vertices of the corresponding basis graph. The 
fifth equality in the equation above is obtained by apply- 
ing to the two 3-cocycles in the fourth row the following 
3-cocycle condition. 

a ([12], [23], [34]) a ([12'], [2'2], [23]) a ([2'2], [23], [34]) 



a([12'],[2'3],[34]) 



a([12'],[2'2],[24]) 



We now explain why this 3-cocycle condition is applica- 
ble. The action of T3 on the basis graph in the first row of 
Eq. (|B10j) followed by an action of T 2 on the resulted ba- 
sis vector in the third row of the equation can be viewed 
(see Fig. [S] and assume the line [2'2] nonexisting for the 
moment) as if the original tetrahedron [1234] is first flat- 
tened to be the triangle [134] and then lifted again to a 
new tetrahedron [12'34]. These two tetrahedra share a 
face, i.e., the triangle [134], so there are seven triangles 
all told. The restriction to the subspace 'H Bf ~ 1 imposes 
a chain rule of the three group elements on each of the 
seven triangles, such that the following identities on the 
group elements hold. 

[2'1] • [12] = [2'3] • [32], 
[2'3] • [32] = [2'4] • [42], 
[2'4] • [42] = [2'1] ■ [12]. 

This enables one to add the line [2'2] and assign [2' 2] = 
[2T] • [12] = [2'3] • [32] = [2'4] • [42]. Hence there are 
now five tetrahedra in Fig. [5] with Bf = 1 on any of the 
ten triangles, justifying the 3-cocycle condition above. 
Furthermore, [2'2] determines [12'], [2'3] and [2'4], the 
summation in the fifth equality in Eq. ^ can be replaced 
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by J2[2'2] /I^L where the factor of 1/\G\ arises because 
[2'2] € G is a new element to be summed over. 




These four identities leads to 
I x (a,b) = 



a{ab,x ,x)a{x ,x,x x ) 



a(a,b, x 1 )a(b,x 1 ,xax 1 )a(x 1 ,xax 1 , xbx 1 ) 

(C2) 

Using ab = ba, We evaluate I x (a, b)/I x (b, a) directly and 
arrive at Eq. (|3"Tj) . 



FIG. 9: The topology of the action of T 2 T 3 . 

We infer from Eq. (0) that T 2 T 3 = 1 on the ground 
states, which together with T 3 T 2 = 1 on H Bf=1 , implies 
that T 2 = Tl and T 3 = t\ on H°. That is, T 2 and T 3 are 
unitary on the ground states. 

An interesting byproduct of this proof is that on the 
subspacc H Bf=1 , 

T 2 T 3 =A V , (Bll) 
where v is the vertex annihilated by the action of T3. 

Appendix C: Ground state operators 

We prove Eq. (f3"Tj) in this appendix. We start with 
the coefficient in Eq. I|28p . which can be expressed by 
I x (a,b)/I x (b,a), where I x (a,b) is 

I x (a, b) = a(a, 6a; -1 , x)a(bx~ 1 , x, ax^ 1 )a(x, ax^ 1 , xbx^ 1 ) 

(CI) 

We use the 3-cocycle condition ([3|) together with 
the normalization condition ^ to rewrite the three 3- 
cocycles above. By 

Sa(a, b, x' 1 , x) = 1 
Sa(b, x^ 1 , x, ax' 1 ) = 1 
5a(x,x~ ,xax~ 1 ,xbx~ 1 ) = l 
5a(x _1 , x, xax~ x ) = 1 



and have 

a(a, x) 



a{bx 1 ,x,ax ) 
a(x, as -1 , xbx^ 1 ) = 



a{ab,x ,x) 
a(b, x~ l , x)a(a, b, x~ r ) 

a(x~ x , x, ax~ 1 )a(b, x~ l , x) 



)i(b,x l ,xax - 1 ) 



1 



a(x l ,xax 1 ,xbx x )a{x,x 1 ,xax 1 ) 
a(x _1 , x, as -1 ) = a(x, xax~ 1 )a(x~ 1 , x, 



Appendix D: Modular Transformations 

Here we show how we construct the modular transfor- 
mation operators S x and T x ■ 



S x : 



^a([2'l],[13],[34])a([2'l],[12],[24])- 1 

^a([2'l],[13],[34])a([2'l],[12],[24])- 1 
2 ^_2' 

ca([2'4'],[4'l],[13])-" 



1 



4' 



a([2'l],[13],[34]) 



a([2'l],[12],[24])a([2'4'],[4'l],[13]) 
T __2' 

'o'i ro'i 



;a ([l'2'],[2'l],[12])- 



1 



4' 



^a([2'l],[13],[34])a([2'l],[12],[24])- 1 

Xa{[2'^%m- l a{[l'2\[2'\Ul2])- 1 

I' 2' . 

a([2'3'],[3'4'],[4'l]) I I \ 



a([l'2'],[2'3'],[3'l]) I I / 

3' 4' 

^a([2'l],[13],[34])a([2'l],[12],[24])- 1 

xa([2'4'],[4'l],[13]r 1 «([l'2'],[2'l],[12])" 1 

xa([2'3'],[3'4'],[4'l])a([l'2'],[2'3'],[3'l])- 1 

1' 2' 

xa([l'2'],[2'3'],[3'4']) J\ J \ 

3' 4' 



2' 



(Dl) 
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T' 



^a([2'l],[13],[34])a([2'l],[12],[24])- 

^a([2'l],[13],[34])a([2'l],[12],[24])- 1 

2 2' 

xa([3"2'],[2'l],[13]) V/* \ 

1 3" 
^a([2'l],[13] ) [34])a([2'l],[12],[24])- 1 
4" 2' . 



a([3"2'],[2'l],[13]) 



■([4"2'],[2'1],[12]) 



3" 



^a([2'l],[13],[34])a([2'l],[12],[24])- 1 
xa([3"2'],[2'l],[13])a([4"2'],[2'l],[12]) 



a([l"3"],[3"2'],[2'l]) 







1" 


3 



2' 



([1"4"],[4"2'],[2'1]) 



^«([2'l],[13],[34])a([2'l],[12],[24])- i 
xa([3"2'] ! [2'l],[13])a([4"2'],[2'l],[12])- 1 
xa([l"3"], [3"2'], [2'1]) a ([l"4"], [4"2'], [2'1])" 



xa([l"3"],[3"4"],[4"2']) 





\ 1 


1—1 


3 



(D2) 



We verify that p xgx 



is indeed a /3 xgx -i represen- 
tation. We check the the multiplication rule for all 
hi, hi £ Z g 



P^hi)^) 

P g {hi,h 2 )p%{hih 2 ) 

P g (hi,h 2 ) 



p^*-\xhix- 


l )p^~' \xh 2 x- 1 ) 




/3 g (x _1 , xhix 


_1 ) /3 g (x _1 ,xh 2 x~ 


l ) 




Pg{h2iX- r ) 




j3 g (x , xhix 


_1 ) P g (x~ 1 ,xh 2 x~ 


') 




Paih^x- 1 ) 




/3 g (x _1 , xhix 


_1 ) f3 g (x~ 1 , xh 2 x~ 






p g (h 2 ,x- r ) 





Pg(x 1 ,xh 1 h 2 x x ) 



p g {hih 2 ,x r ) 



p'f \xh 1 h 2 x- 1 ) 



ri 9 {hi,x)r) a {h 2 ,x) ,^ 9I -i, , , _ x , 

= r,°(hih 2 ,x) P'faMpS (^ih 2 x ) 

=/3 X gx-i{xh 1 x~ 1 ,xh 2 x~ 1 )pZ 9X (xhifox^ 1 ) 



(E2) 



where the last equality is obtained by using the following 
relation: 



rj 9 (hi,x)rj 9 (h 2 ,x) P xgx -i(xhiX 1 ,xh 2 x 1 ) 



rj 9 (hih 2 , x) 



P g (hi,h 2 ) 



(E3) 



which can be verified by directly applying the 
twisted 2-cocycle conditions 8f3 g = 1 in 



ccssively to the triples (hi, h 2 , x^ 1 ), (hi,x~ ,xh 2 x~~ ), 
(x^ 1 , xhix , xh 2 x~ 1% ) and make the corresponding sub- 
stitutions. 

An immediate consequence of the above isomorphism 
is the relation between the projective characters, 



(xhx 



r, 9 (h,x)x 9 Jh), 



(E4) 



which is the very Eq. (|4T))) . From the relation, we obtain 
the following proposition. 



Appendix E: Solutions for 5 and T matrices 

1. Projective Characters. 

The ccntralizer subgroups Z g are isomorphic for all el- 
ements g of the conjugacy class C A . Therefore the corre- 
sponding projective representations are also isomorphic. 
Given the representation ]fi for a fixed g £ C A , we con- 
struct p X9 ' x 1 as follows. For all x £ G, the elements 
xhx -1 runs over all elements in xZ g x~ x while h runs 
over all elements in Z g . We can define a projective repre- 
sentation p X9X of xZ g x~ l from a given representation 
Pi of Z g , by 

^(xhx- 1 ) = ^^frP ^Cfe) (El) 



Proposition 1 If h £ Z g is not (3 g ~regular, X 9 a{h) = 0. 

Proof. The proof is straightforward. If h £ Z g is not /3 g - 
regular, there must exist k £ Z g> h, such that f3 g (h,k) ^ 
{3 g (k, h). By setting x — fc -1 in Eq. (|E4|) we have 



2. Some Proofs 

We prove the equalities in Eq. in particular the 

last one as follows. We first write down the action of 
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T\A, fi) explicitly. 

T\A,fx) =T e P°\A,ri=T e \A :f i) 



r(Z A ,P aA ) 



= W\ ^ 



g£C A ,h£Z g v=l 

g (-rr„Y„-i u-i „ j.n„y„-ii,-i i, 



xa(h,g- 1 h- 1 ,g)a(g,g~ 1 h- 1 ,hy 1 
xa{g,h,g- 1 h~ 1 )a(h,g,g- 1 h- 1 )- 1 



xa(g,hg 1 ,g) X i{g~ 1 h~ 1 )\A,v), 

where the inverse transformation in Eq. (|66p is used. 
Note that g must belong to precisely one conjugacy class, 
so in the equation above the classes C A and C B are iden- 
tified, which kills the summation over C B . Knowing 
this, we rewrite the action of T in the equation above 
by grouping the 3-cocycles into the twisted 2-cocycles. 

r\A,(j) 

=T^ E E ^g- l h)xl{h)xl{g-^h) 

1 1 g ec A "=i 
hez g 



(E5) 



P g (h,g^h^)(3 g {g^h,gh) 
P g (g-ih-\h)p g (gh,g-ih) 1 ' ; 

r(Z A ,P gA ) 

=Tn E E ^g-^ximiig-^h) 

11 g€C A "=1 

=Tn E E ^(5,5-^)^(^)^(5-^)1^^) 

11 gec A "=1 

where in the last step use is made of Eq. P§| to absorb 
the two rj 9 terms into the projective characters x% an d 

A. I/" 

Because Eq. says that jo£(<?) is a multiple of the 
identity matrix, we have 

(3 g {g,g- l h)x%{h) = tr [A^jT^JW] 
= tr[^ (5)^(5-^)] 



tr[ 



x£ (5 A ) , 
dim,, 



1^(5-^)] 



x£ (5 A ) : 
dim,, 



Clearly, = ^2 g -i h , together with the orthogonality 
condition in Eq. (|50[) . the summation evaluates to 



Y gA (q A ) 



(E6) 



confirming Eq. 

In the sequel, we derive the 5-matrix in Eq. (JS3J). We 
first act the S operator on a generic eigenvector 1 73, v) of 
the T operator. 

S\B,v) =S e P°\B,is) =S e \B,is) 
1 



1^1 g'EC B ,h'£Z 



E X 9 J(h') 



xa(5'- 1 /i'- 1 ! 5',/i')a(5'- 1 /i'- 1 ! /i / ,5 / )- 1 



g'rMg'^g'-'h'-^h')- 1 



xaih'-^g'- 1 »> 
xar 1 ^ 1 ,,',/- 1 )^^ 1 ,?- 1 ^ 1 ,^ 1 
xaig^g'-'h'- 1 ,^'- 1 ^'}. 



Notice that in the equation above, the first two 
and the fourth a terms define a twisted 2-cocycle 
0g' (5' -1 ^' -1 ? h')~ l the sixth a equals 1 as it is normal- 
ized, and by the 3-cocycle condition the fifth and seventh 
a terms are equal to a(/i' _1 , g', g'~ 1 )a(g' , g'~ 1 h'~ 1 , g')^ 1 , 
which together with the third a term, define another 
twisted 2-cocycle As such, the above 

equation becomes 



S\B,v) 



E X g u'(h') 



g'GC B ,h'eZ Q 



x/3 g '(g 



h')- 1 f}Ah'- 1 ,s'- 1 r 1 \h'- 1 ,g') 



xZ'(h<-i)\h'-\g'), 



where the second equality and third equality are respec- 
tively the results of the following two identities. 

pi{h')p({h'- r ) = P g -{h'-\h')t 



«• xt(h') = x 9 l/ (h>-i)p 3 ,(h'-\h% 

which is due to the unitarity of the projective represen- 
tation p. 



P g ,{h>-\h')P g ,{g'-\e) 
P g ,(g'-^h>-\h>)P g ,{g>-\h>-^) 

f}g'(h'-\h') 
p g ,{g>-ih'-\hi)p g .{h>-\gi-i) 



= 1, 



which is a consequence of the (twisted) 2-cocycle con- 
dition, the normalization of /3 g ', and the fact that 
/Ms' -1 ) fr' -1 ) = Pg'ih'- 1 ,^' 1 ) because h'- 1 must be 
/3 g '-regular otherwise x 9 {h'^ 1 ) vanishes. 
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As such, the S-matrix reads where (g,h\h' 1 ,g') — <5;i'-i. g <V./i is understood. This 

proves Eq. (|6"3"1) . 

S(Ap,)(Bu) 



{A,n\S\B,v) 

I^E E W(Vx S J(h>-i){g,h\h'-\g') 



\G. 

'gec A 9 'ec B 
hez g h'ez g , 

1 1 g ec A hec B 
h£Z Q g&z h 
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